Sflfi  riU  ^ 


THE  BOM  CORPORATION 


A-FOfcsR-TH-  8  3-0326 


7915  Jones  Branch  Drive.  McLean.  Virginia  22102-3396  •  (703)  848-5000  ♦  Telex;  901103  BDM  MCLN 


A  New  Perspective  on 
Rutting  in  Flexible  Pavements 


Final  Technical  Report 


PREPARED  FOR  AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH,  BOLLING  AIR 
FORCE  BASE,  WASHINGTON,  D.C.  20332-6448 


DTIC. 

ELECTEg? 
MAR  2  3  198^Ju 

1  c*  i 


FEBRUARY  24,1989 


BDM/MCL-89-0108-TR 


A<  -rc7,:d  rubllc  riltKKo; 

’  89 

3  1 

bdtb 

THE  BDM  CORPORATION 


791 5  Jones  Branch  Drive.  McLean.  Virginia  22102-3396  •  (703)  848-5000 


A  NEW  PERSPECTIVE  ON  RUTTING  IN 
FLEXIBLE  PAVEMENTS 
FINAL  TECHNICAL  REPORT 
February  24,  1989 


BDM/MCL-89-0108-TR 


Prepared  for  Air  Force  Office  of  Scientific  Research,  Bolling  Air  Force 
Base,  Washington,  OC  20332-6448. 


■KIP  1  I  II II  MM 

f:nsMrwnr*ir.T7r.T:i 


REPORT  DOCUMENTATION  PAGE 

tj  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

lb  RESTRICTIVE  MARKINGS 

2a  SECURITY  CLASSIFICATION  AUTHORITY 

3  DISTRIBUTION/  AVAILABILITY  OF  REPORT 

Approved  for  Public  Release; 

Distribution  Unlimited 

2b  DECLASSIFICATION  -  DOWNGRADING  SCHEDULE 

4  PERFORMING  ORGANIZATION  REPORT  NUMSER(S) 

5  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

TR-  8  9  -  0  326 

6a  NAME  OF  PERFORMING  ORGANIZATION 

The  8DM  Corporation 

6b  OEFICE  SYMBOL 
(if  opplicsbto) 

7a.  NAME  OF  MONITORING  ORGANIZATION  ; 

AFOSR/NA  | 

6c  ADDRESS  (G'fy.  Sf4f*.  ond  ZIP  Cod*) 

7915  Jones  Branch  Drive 

McLean,  VA  22102 

7b  AOORESS  (Oty,  Stott,  tnd  ZIP  Code) 

Building  410 

Bolling  AFB,  DC  20332-6448 

8a.  NAME  OF  FUNDING  /  SPONSORING 
ORGANIZATION 

AFOSR 

8b  OFFICE  SYMBOL 
(If  spptKsbte) 

NA 

9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F4962088-C-0019 

Building  410 

Bolling  AFB,  DC  20332-6448 


1 1  title  (Includ*  Srtunty  aoaifkttion) 

A  New  Perspective  On  Rutting  In  Flexible  Pavements 


12  PERSONAL  AUTHOR(S) 

Dr.  K.  0.  Harrop-Williams 


PROJECT 

TASK 

NO 

NO 

2302 

Cl 

13a.  type  Of  REPORT 

Final 


16  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVEREO 
PROM  12/15/B7TO 


14  DATE  OE  REPORT  (Yor.  Mont*.  Oay) 


YJt  w.w&rjM  Kray,  MB  n  T^rrm  H 


COSATI  COOES 


GROUP 


H.  SUBJECT  terms  {CormiHgi  on  rovorto  if  nornuty  ond  idtntity  6> 

evrs,.  3 

cot/riJoiy^,  +  ft.  fHfr-rz  i E 


SUB-GROUP 


>  on  reverse  1/  necessary  and  identify  by  btock  number) 

e/v/rs,  3  7R& 


£v  /f-4 


T9  ABSTRACT  {Continue  on  reverse  if  necesury  end  identify  by  btock  number) 

This  report  presents  a  new  approach  to  the  prediction  of  rutting  in  flexible 
pavements.  It  recognizes  the  prominent  role  played  by  the  granular  layer  in  the  stress 
distribution  and  strain  accumulation  in  flexible  pavements  and  diverts  from  the 
conventional  continuum  approach  to  modeling  stresses  in  the  layer.  Fundamentally,  the 
conventional  requirement  for  the  existence  of  the  second  derivatives  of  strains  in  the 
granular  layer  is  replaced  by  a  stress  continuity  equation  that  provides  for  stress 
transfer  between  particles  through  these  contacts.  This  is  developed  assuming  that 
deformation  is  due  primarily  to  particle  movements  and  is  shown  to  be  based  on  a  nonlinear 
stress-strain  relationship  unique  for  granular  material.  In  order  to  evaluate  the  rutting 
in  a  flexible  pavement,  the  particulate  theory  of  stress  transfer  is  extended  to  predict 
stresses  in  multilayered  systems  consisting  of  granular  and  elastic  layers.  Knowledge 
of  the  nature  of  stress  transfer  and  stress-strain  response  is  next  incorporated  into 
a  theory  for  strain  accumulation  with  repetitive  loading. 


20  DISTRIBUTION  /AVAILABILITY  OP  ABSTRACT 

Qunclassipieo/unlimiteo  □  same  as  ret 


00  FORM  1473, 34  MAR 


121  ABSTRACT  SECURITY  CLASSIFICATION 

UNCLASSIFIED  


22c.  OFFICE  SYMBOL 

AFOSR/NA 


83  APR  edition  may  Qe  used  until  exhausted 
All  other  editions  are  obsolete 


JRITY  CLASSIFICATION  Of  THi$  PAGE 

UNCLASSIFIED 


1 


preface 


This  report  was  prepared  by  The  BDM  Corporation.  7915  Jones  Branch 
Drive,  McLean.  Virginia  22102  under  Contract  Number  F4962088-C-0019  for  the 
Air  Force  Office  of  Scientific  Research,  Bolling  Air  Force  Base, 
Washington,  D.  C.  20332-6448. 

This  report  summarizes  work  done  between  December  15,  1987,  and 
February  15,  1989.  AFOSR  Program  Manager  was  Major  Steven  Boyce,  USAF. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  (PA)  and 
is  releasable  to  the  National  Technical  Information  Service  ( NT I S ) .  At 
NT  I S ,  it  will  be  available  to  the  general  public,  including  foreign 
nati onal s . 

This  technical  report  has  been  reviewed  and  is  approved  for 
publ i cati on . 


Steven  Boyce,  Major  (USAF) 

Program  Manager 

Directorate  of  Aerospace  Sciences 


A  r  ■■ 
N ,  . 


on  For 
'•i  I 


~vT 


;  l 


table  of  contents 


Section  Title  Page 


PREFACE .  iii 

TABLE  OF  CONTENTS .  iv 

LIST  OF  FIGURES .  vi 

I  INTRODUCTION .  1 

A.  OBJECTIVES .  1 

B.  BACKGROUND .  2 

C.  SCOPE .  5 

d!  ORGANIZATION  OF’TECHNiCAL  REPORT]  1  1  6 

II  STRESS  TRANSFER  IN  GRANULAR  MATERIALS  UNDER  INCLINED 

LOADS .  7 

A.  INTRODUCTION .  7 

B.  TRADITIONAL  APPROACH  TO  STRESS  DISTRIBUTION  IN 

GRANULAR  MEDIA .  8 

C.  NEW  APPROACH  TO  STRESS  DISTRIBUTION  IN  GRANULAR 

MEDIA .  9 

1.  Two  Dimensions .  13 

2.  Three  Dimensions .  14 

D.  SOLUTION  FOR  STRESSES  UNDER  INCLINED  LOADS .  15 

1.  Two  Dimensions .  15 

2.  Three  Dimensions .  17 

III  STRAIN  ACCUMULATION  IN  GRANULAR  MATERIALS  UNDER  STATIC 

AND  REPEATED  LOADS .  18 

A.  INTRODUCTION .  18 

B.  TRADITIONAL  APPROACH  TO  STRESS - STRAI N  MODELING  OF 

GRANULAR  MATERIAL .  18 

C.  NEW  APPROACH  TO  STRESS-STRAIN  MODELING  OF  GRANULAR 

MATERIAL .  20 

1.  General  Theory .  20 

2.  Application  To  One  Dimensional  And  Isotropic 

Loading  Conditions .  21 

3.  Application  To  Triaxial  Loading  Conditions..  22 

D.  PERMANENT  AXIAL  STRAIN  ACCUMULATION  IN  GRANULAR 

MATERIALS .  24 

1.  Static  Loading .  24 

2.  Repeated  Loading .  25 

E.  PERMANENT  RADIAL  STRAIN  ACCUMULATION  IN  GRANULAR 

MATERIALS .  28 


iv 


TABLE  OF  CONTENTS 
(CONCLUDED) 


Section  Title  Page 

IV  VALIDATION  OF  PARTICULATE  THEORY  FOR  GRANULAR  MATERIAL.  30 

A.  INTRODUCTION .  30 

B.  VALIDATION  OF  PARTICULATE  THEORY  OF  STRESS 

TRANSFER .  30 

C.  VALIDATION  OF  PARTICULATE  iriEORY  OF  STRESS- STRAI N 

RESPONSE .  35 

D.  VALIDATION  OF  THE  PARTICULATE  THEORY  OF  STRAIN 

ACCUMULATION  UNDER  REPEATED  LOADING .  39 

V  EVALUATION  OF  MATERIAL  PARAMETERS  IN  PARTICULATE  THEORY  43 

A.  INTRODUCTION .  43 

B.  MATERIAL  CONSTANTS  IN  PARTICULATE  THEORY  OF 

STRESS  TRANSFER .  43 

C.  MATERIAL  CONSTANTS  IN  PARTICULATE  THEORY  OF 

STRESS-STRAIN  RESPONSE .  44 

D.  MATERIAL  CONSTANTS  IN  PARTICULATE  THEORY  OF 

STRAIN  ACCUMULATION .  47 

VI  RUTTING  MODEL  FOR  MULTILAYERED  PAVEMENTS  WITH 

PARTICULATE  MATERIAL .  49 

A.  INTRODUCTION .  49 

B.  STRESS  TRANSFER  IN  MULTILAYERED  MEDIA .  50 

C.  RUTTING  PREDICTION  MODEL .  53 

VII  CONCLUSIONS  AND  RECOMMENDATIONS .  56 

A.  CONCLUSIONS .  56 

B.  RECOMMENDATIONS .  57 

LIST  OF  REFERENCES .  59 

APPENDIX 

A  A  GENERAL  STRESS-STRAIN  MODEL  FOR  GRANULAR  SOILS .  62 


v 


LIST  OF  FIGURES 


Figure  Title  Page 

1  Adjacent  Soil  Elements .  12 

2  Vertical  Stress  Below  The  Center  Of  A  Uni formly- Loaded 

Ci  rcul  ar  Area .  32 

3  Location  Of  A  Point  P  From  Circular  Loaded  Area  Of 

Radius  A .  33 

4  Stress-Strain  Behavior  For  Isotropic  Loading .  36 

5  Stress-Strain  Behavior  For  Triaxial  Compression .  37 

6  Relationship  Of  Volumetric  And  Axial  Strains .  38 

7  Log-Log  Plot  Of  Permanent  Axial  Strain  Versus  Number 

Of  Load  Cycl  es .  40 

8  Semi -Log  Plot  Of  Permanent  Axial  Strain  Versus  Number 

Of  Load  Cycles .  41 

9  Relation  of  E0.  £u.  And  £e  To  Confining  Pressure .  46 


vi 


SECTION  I 
INTRODUCTION 


A.  OBJECTIVES 

The  objective  of  this  research  is  to  develop  basic  theories  that 
can  be  used  to  predict  rutting  in  flexible  pavements  with  granular 
layers.  The  approach  taken  in  this  research  is  to  treat  the  granular 
material  as  particulate  in  nature  ana  in  this  way  provide  a  more 
rational  model  than  conventional  theory  that  treats  the  particulate 
material  as  a  continuum.  It  takes  a  premise  that  the  major  reason  for 
the  inability  of  conventional  elastic  or  elasto-plastic  theories  to 
model  granular  soil  behavior  is  the  fact  that  these  soils  consist  of 
randomly-arranged,  i rregularly-shaped  discrete  particles  that  are  free 
to  displace  relative  to  each  other.  It  assumes  that  the  mechanism 
responsible  for  the  deformation  of  the  material  subject  to  induced 
loading  is  not  caused  so  much  by  deformation  of  individual  particles  as 
by  relative  movements  of  the  more  mobile  particles  forming  a  denser 
state . 

In  linear  and  nonlinear  continuum  theory  an  equation  ensuring  the 
existence  of  the  second  derivatives  of  strains  throughout  the  media  is 
introduced.  This  equation,  through  constitutive  relations,  is  combined 
with  the  equilibrium  equations  for  stresses  to  provide  problem  solution. 
One  inaccuracy  in  using  this  theory  to  model  the  behavior  of  granular 
material  is  the  inappropriateness  of  this  strain  requirement.  This 
requirement  inevitably  results  in  the  prediction  of  horizontal  tensile 
stresses  when  the  granular  media  is  subjected  to  inclined  loads.  It 
will  also  predict  stresses  if  the  granular  media  is  subjected  to  purely 
horizontal  surface  loads,  an  unrealistic  condition  since  cohesionless 
material  will  provide  no  resistance  to  these  loads.  This  is  the  reason 
why  the  continuing  development  of  sophisticated  nonlinear  empirical 
equations  to  model  the  experimentally-observed  stress- strai n  relations 
of  the  soil  does  not  greatly  improve  the  continuum  mechanics  prediction. 
Khedr  (1985)  observed  that  even  finite  element  methods  using  elasto- 
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plastic  theory  predicts  unrealistic  radial  pressure  because  they 
characterize  the  granular  layer  as  a  continuum. 

In  the  particulate  theory  developed  here  the  requirement  for  the 
existence  of  the  second  derivatives  of  strains  in  the  granular  media  is 
replaced  by  a  stress  continuity  equation  that  guarantees  stress  transfer 
between  particles  through  their  contacts.  This  continuity  equation  is 
derived  from  two  points  of  view,  namely:  (1)  from  the  assumption  that 
de forma t i on  is  caused  by  particle  movements  rather  than  particle 
compression;  and  (2)  that  in  the  definition  of  stress  and  strain,  any 
represemati ve  element  of  soil  must  be  composed  of  particles  and  voids. 
In  this  way.  the  elemental  volume  of  soil  cannot  be  made  to  approach 
zero  but  must,  instead,  approach  a  finite  minimal  volume  with  enough 
particles  such  that  the  particle  movements  produce  strain.  The 
combination  of  this  continuity  equation  with  the  equilibrium  equation 
allows  the  determination  of  all  the  components  of  stress.  In  addition, 
the  particulate  theory  results  in  a  stress-strain  response  model  unique 
for  granular  soils. 

The  specific  objectives  of  this  research  are  to: 

(1)  Develop  a  particulate  theory  to  predict  stress  transfer 
through  granular  material  subjected  to  inclined  loads. 

(2)  Extend  the  particulate  theory  to  predict  the  stress- strai n 
response  under  static  loading; 

(3)  Use  the  particulate  theory  to  predict  strain  accumulation 
under  repeated  loading;  and 

(4)  Combine  the  results  of  objectives  (1),  (2).  and  (3)  to 
predict  rutting  in  mul ti 1 ayered ,  flexible  pavements  with  a  granular 
1 ayer . 

8.  BACKGROUND 

Flexible  pavement  design  requires  the  ability  to  predict  pavement 
performance.  One  of  the  major  indicators  of  pavement  performance  is  the 
distress  caused  by  rutting.  This  rutting  is  the  accumulation  of 
permanent  settl ement  with  traffic  application,  and  for  flexible 
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pavements  it  occurs  in  all  inelastic  components  of  the  pavement. 
Rutting  will  eventually  result  in  a  reduction  of  pavement  serviceability 
due  to  less  in  riding  comfort  and  may  also  lead  to  hydroplaning  and 
icing  due  to  the  collection  of  the  water  in  the  deformati  ons .  As  a 
result,  any  rational  method  of  predicting  rutting  in  flexible  pavement 
must  oe  aoie  to  model  both  the  stress  transfer  mechanism  and  the  stress- 
strain  characteristics  of  each  layer.  It  must  also  be  able  to  combine 
the  behavior  of  each  layer  -;n  the  way  that  it  contributes  to  the  overall 
performance  of  the  pavement. 

It  is  common  practice  in  pavement  design  to  assume  stress 
distribution  using  the  multiple  elastic  layers  theory  because  of  the 
relative  ease  at  whic.i  solutions  may  be  obtained  (Monismith  and  Finn. 
1977;  Monismith.  Finn,  and  Epps.  1986).  The  derived  stresses  which  are 
based  on  linear  stress-strai n  relationships  that  acknowledge  no  strain 
accumulation  are  next  used  with  empirical  or  el asto-plastic  nonlinear 
stress-strain  relationships  to  predict  strain  and  strain  accumulation 
with  repeated  loading.  However,  experimental  data  show  that  the 
response  of  soil,  and  especially  granular  materials,  depends  strongly  on 
the  state  of  stress.  Therefore,  it  is  essential  to  know  the  correct 
stress  conditions  in  the  flexible  pavements  before  strain  accumulation 
can  be  rationally  predicted.  Most  of  the  research  in  the  last  decade 
has  concentrated  on  developing  stress  prediction  that  includes  the 
nonlinearity  of  the  stress - stra i n  relationship  of  granular  material. 
This  usually  takes  the  form  of  empirically  relating  the  resilient 
modulus,  defined  as  the  ratio  of  the  repeated  deviator  stress  to  the 
recoverable  strain,  to  the  sum  of  the  principal  stresses.  Besides 
requiring  complex  finite  element  solutions,  this  model  has  serious 
limitations.  More  complex  and  soph i s t i cated  models  giving  better 
descriptions  of  resilient  response  do  not  greatly  improve  the  prediction 
(llzan,  1985;  Brown  and  Pappin.  1981). 

It  is  very  difficult  to  accurately  predict  rutting  in  flexible 
pavement  because  the  granular  layer  is  the  most  important  load-carrying 
component  of  the  pavement.  To  adequately  design  flexible  pavements,  an 
accurate  understanding  of  the  behavior  of  granular  material  is 
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necessary.  This  is  more  essential  now  than  ever  since  existing 
pavements  and  future  pavements  are  expected  to  accommodate  heavier  and 
heavier  loads.  This  material  consists  of  discrete  particles,  and  it 
deforms  as  the  result  of  particle  movements  rather  than  particle 
deformation.  Experimental  evidence  abounds  (Haggarty,  1963;  Morgan  and 
Gerrard.  1971)  f  or  the  inaD.lity  of  classical  elastic  or  inelastic 
continuum  mechanicc  to  predict  granular  materials  behavior  under  load. 
Tne  main  reason  for  this  is  tnat  even  the  most  sophisticated  of  these 
solutions  requires  the  second  derivatives  of  strains  to  exist  at  all 
points  in  the  granular  media.  This  requirement  is  unrealistic  as  voids 
always  exist  in  the  material,  and  tne  deformations  are  the  result  of 
discrete  particle  displacements. 

An  alternative  approacn  to  stress  distribution  in  granular 
materials  can  also  be  found  in  the  literature  (Golden,  1984,  1986;  Harr, 
1977;  Hill  and  Harr.  1982;  Endley  and  Peyrot,  1977;  Chikwendu  and 
Alimba,  1979;  Sergeev.  1969).  This  approach  considers  the  discrete 
nature  of  the  soil  and  assumes  that  when  a  normal  point  force  is  applied 
on  the  surface  of  an  infinite  half-space,  the  influence  of  this  force 
travels  from  particle  to  particle  in  a  fashion  analogous  to  a  random 
walk  or  continuous  Markov  process.  The  resulting  equation  is  the 
diffusion  equation  for  vertical  stresses.  This  approach  ignores  the 
requirement  for  the  existence  :f  the  second  derivatives  of  strains  and 
requires  knowledge  of  the  diffusion  coefficient,  a  material  property. 
The  results  of  this  approach  show  that  it  provides  a  better  qualitative 
fit  to  observed  stress  distributions  in  granular  media  and  can  be  used 
to  predict  the  stresses  in  layered  media  if  the  diffusion  coefficients 
of  the  layers  are  known  (Harr,  1977;  Golden.  1984).  The  limitations  of 
this  theory  at  this  point  are  the  description  of  the  diffusion 
coefficient  and  the  role  the  constitutive  nature  of  the  material  plays 
in  the  stress  diffusion. 

Even  from  a  conceptual  viewpoint,  any  particulate  theory  is  an 
improvement  over  continuum  models.  As  granular  soils  are  particulate  in 
nature,  stresses  are  transmitted  at  particle  contacts,  and  strains  are 
the  result  of  particle  movements  rather  than  particle  compression.  The 
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vertical  stress  distribution  predicted  by  the  stochastic  theory  for  a 
point  load  on  the  surface  agrees  with  experimental  results  which  show  a 
more  bell-shaped  stress  distribution  than  the  elastic  prediction  (Harr, 
1977).  Any  particulate  theory  based  on  stress  diffusion  from  a  source 
should  give  a  bell-shaped  distribution  as  this  is  the  solution  of  the 
diffusion  equation  for  a  point  source.  Further,  with  the  requirement 
for  the  specification  of  a  material  property  in  the  form  of  the 
diffusion  coefficient,  the  stochastic  theory  shows  that  stress  transfer 
is  dependent  on  particle  sizes,  shapes,  packing,  load  history,  etc. 
This  is  in  opposition  to  the  elastic  theory  which  predicts  the  same 
vertical  stress  at  a  given  point  in  sand  as  it  does  in  steel. 

C.  SCOPE 

This  report  presents  a  new  approach  to  the  prediction  of  rutting 
in  flexible  pavements.  It  recognizes  the  prominent  role  played  by  the 
granular  layer  in  the  stress  distribution  and  strain  accumulation  in 
flexible  pavements  and  diverts  from  the  conventional  continuum  approach 
to  modeling  stresses  in  the  layer.  The  methodology  of  the  new  approach 
is  as  follows: 

(1)  The  requirement  for  the  existence  of  the  second  derivatives 
of  strains  in  the  granular  layer  conventionally  used  in  the 
determination  of  stress  is  relaxed.  It  is  replaced  by  a  stress 
continuity  equation  that  provides  for  stress  transfer  between  particles 
through  these  contacts.  This  is  developed  assuming  that  deformation  is 
due  primarily  to  particle  movements: 

(2)  It  is  shown  that  the  development  of  the  stress  continuity 
equation  is  based  on  a  nonlinear  stress- strai n  relationship  unique  for 
granular  material.  This  relationship  is  derived  by  the  particulate 
approach  and  reduces  to  the  hyperbolic  model  for  soils  that  do  not 
exhibit  a  distinct  peak  stress; 

(3)  Knowledge  of  the  nature  of  stress  transfer  and  stress-strain 
response  in  granular  media  is  next  incorporated  into  a  theory  for  strain 
accumulation  with  repetitive  loading:  and 
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(4)  In  order  to  evaluate  the  rutting  in  a  flexible  pavement,  the 
particulate  theory  of  stress  transfer  is  extended  to  predict  stresses  in 
multi-layered  systems  consisting  of  granular  and  elastic  layers.  In 
this  way,  more  accurate  predictions  of  the  stresses  in  the  layers  will 
lead  to  more  accurate  predictions  of  rutting. 

D.  ORGANIZATION  OF  TECHNICAL  REPORT 

This  technical  report  is  divided  into  eight  sections  and  an 
appendix.  Section  I  is  the  introduction  which  outlines  the  objectives, 
background  and  scope  of  the  research.  Section  II  investigates  the 
stress  transfer  mechanism  in  granular  media.  After  a  look  at  the 
traditional  methods,  a  new  particulate  approach  to  stress  transfer  is 
introduced.  Section  III  presents  a  particulate  theory  for  stress-strain 
response  in  granular  materials.  It  addresses  both  the  stress-strain 
response  due  to  static  loading  and  permanent  strain  accumulation  under 
repeated  loading.  Section  IV  presents  the  validation  of  the  theories 
developed  in  Sections  II  and  III  using  information  found  in  the 
literature.  Methods  of  determining  the  experimental  constants 
identified  in  the  theories  of  Sections  II  and  III  are  presented  in 
Section  V.  Section  VI  presents  a  method  for  predicting  rutting  in 
flexible  pavements  with  granular  layers.  Section  VII  presents  the 
Conclusion  and  Recommendations,  and  Section  VIII  is  the  List  of 
References.  Finally,  a  paper  published  from  this  work  is  presented  in 
an  appendix. 
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SECTION  II 

STRESS  TRANSFER  IN  GRANULAR  MATERIALS  UNDER  INCLINED  LOADS 
A.  INTRODUCTION 

Granular  soils  generally  consist  of  randomly  arranged, 
i  rregul  arly- shaped ,  discrete  particles  that  are  free  to  displace 
relative  to  each  other.  The  deformation  of  this  material  under  load  is 
not  caused  so  much  by  deformation  of  individual  particles  as  by  relative 
movement  of  the  more  mobile  particles  forming  a  denser  state.  As  a 
result,  experiments  have  shown  that  deformation  predicted  by  the  theory 
of  elasticity  is  incorrect  in  magnitude  and  distribution  (Turnbull. 
Maxwell,  and  Ahlvin.  1961;  Morgan  and  Gerrard.  1981).  Really  obvious 
di screpanci es  between  experiments  and  elastic  predictions  are  that 
horizontal  tensile  stresses  are  predicted  under  inclined  loads  and 
stresses  are  predicted  under  purely  horizontal  surface  loads.  In 
actuality,  tensile  stresses  cannot  exist  in  cohesionless  materials  and 
cohesionless  materials  will  provide  no  resistance  to  purely  horizontal 
surface  loads. 

The  main  reason  for  the  elastic  theory's  prediction  of  tensile 
stresses  in  granular  materials  under  inclined  loads  and  the  prediction 
of  stress  transmission  in  granular  material  under  purely  horizontal 
surface  loading  is  the  requirement  for  the  second  derivatives  of  strains 
to  exist  at  all  points  in  the  material.  In  elastic  predictions,  using 
linear  stress - strai n  relationships,  this  condition  is  combined  with  the 
equilibrium  equations  to  give  the  stresses  in  the  media.  Traditional 
methods  of  improving  the  predictions  of  stresses  in  granular  methods 
introduce  more  realistic  nonlinear  stress-strain  behavior.  However,  the 
requirement  for  the  existence  of  the  second  derivatives  of  strains  at 
all  points  in  the  media  is  still  imposed.  It  seems  unlikely  that  these 
methods  of  improving  the  prediction  of  stresses  in  granular  material  can 
avoid  the  condition  of  predicting  tensile  horizontal  stresses  under 
inclined  load  or  the  transmission  of  stresses  under  purely  horizontal 
loads . 


An  alternative  approach  to  stress  distribution  in  granular  media 
that  does  not  require  the  existence  of  the  second  derivatives  of  strains 
at  all  points  but  on  stress  transfer  through  particle  contacts  can  be 
found  in  References  10  and  13.  Their  approach  assumes  that  the 
influence  of  a  surface  force  travels  from  particle  to  particle  in  a 
fashion  analogous  to  a  continuous  Markov  process.  The  resulting 
equation  is  the  diffusion  equation  for  vertical  stresses.  However, 
there  is  no  constitutive  relationship  identified  in  the  stochastic 
approach  (Reference  10). 

In  this  section,  the  stress  transfer  mechanism  for  an  inclined 
load  on  a  granular  material  is  developed  by  recognizing  that  particles 
are  of  finite  size  and  that  stress  and  strain  cannot  be  defined  at  a 
point  in  granular  material  but  only  with  respect  to  some  finite  volume. 
It  is  shown  that  this  approach,  although  entirely  deterministic,  leads 
to  the  diffusion  equation  of  the  stochastic  approach.  More  importantly, 
it  shows  that  there  is  a  built-in  nonlinear  stress-strain  relationship. 
The  approach  avoids  the  requirement  for  the  existence  of  the  second 
derivatives  of  strains  and  avoids  the  prediction  of  tensile  stresses  in 
granular  soils  and  the  prediction  of  stresses  under  purely  horizontal 
surface  loads. 

B.  TRADITIONAL  APPROACH  TO  STRESS  DISTRIBUTION  IN  GRANULAR  MEDIA 


The  solution  of  stresses  in  any  body  is  obtained  by  solving  the 
equations  of  equilibrium.  In  two  dimensions,  neglecting  the  weight  of 
the  soil,  these  equilibrium  equations  are 


and 


d Ox  d  Tzx 

3x  +  dz 

doz  dxxz 
'Z  +  dx 


txz  “  xzx 


=  0 

=  0 


(1) 

(2) 


(3) 
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where  az  and  ox  are  the  normal  stresses  on  the  element  in  the  vertical 
and  horizontal  directions,  respect i vely ;  and  tx2  and  tzx  are  the 
vertical  and  horizontal  shear  stresses.  Since  these  are  only  three 
equations  for  the  four  unknown  stresses,  a  fourth  equation  is  necessary 
for  tne  solution  of  the  stresses.  Traditionally,  this  fourth  equation 
is  supplied  by  assuming  the  existence  of  the  second  derivatives  of 
strains  at  all  points  in  the  granular  media.  This  requirement  results 
in  the  compatibility  equation  (Timoshenko,  1951) 

32£x  d2£z  _  a^y 

3x2  +  az2  "  3x3z  (4) 

where  ex  and  e2  are  axial  strains  and  y  is  the  shear  strain.  The 
stresses  and  strains  are  then  related  by  an  appropriate  constitutive 
relationship.  For  the  special  case  of  linear  elastic  assumption, 
closed-form  solutions  for  the  stresses  are  available.  These,  however, 
fail  to  appropriately  model  the  observed  conditions  in  granular  media 
(Morgan  and  Gerrard.  1981).  With  the  use  of  empirical  stress-strain 
models,  such  as  the  hyperbolic  model,  intense  numerical  modeling  is 
necessary.  The  complexity  of  solution  and  the  assumption  of  the 
existence  of  the  second  derivatives  of  strains  at  all  points  in  the 
media  are  the  limitations  of  the  approach. 

C.  NEW  APPROACH  TO  STRESS  DISTRIBUTION  IN  GRANULAR  MEDIA 

The  new  approach  considers  any  representati ve  element  of  soil  of 
volume  dxdydz  tc  be  composed  of  particles  and  must  satisfy  the 
equilibrium  equations.  However,  rather  than  assuming  the  existence  of 
the  second  derivatives  of  strains  to  exist  at  all  points  in  the  granular 
media,  the  particulate  nature  of  soil  deformation  is  examined  to  develop 
the  additional  equations  necessary  for  solution  of  the  stresses. 

Since  any  representati  ve  element  of  soil  must  be  composed  of 
particles  and  voids,  its  volume  dxdydz  cannot  be  made  to  approach  zero 
but  must  approach  a  finite  limiting  volume,  say  ijh,  with  enough 


Q 


particles  so  that  the  relative  movements  of  particles  in  ijh  caused  by 
forces  on  these  particles  produce  strain.  If  w  is  the  average  change  in 
vertical  displacement  of  particles  in  the  element,  then  the  vertical 
strain  in  the  element  is 

ez  =  1  i m  w/dz  =  w/h 

dz^h  (5) 

Also,  if  F  is  the  vertical  component  of  force  on  a  horizontal  plane 
through  the  element,  then  the  vertical  stress  is 

cr2  =  lim  F/(dxdy)  =  F/(ij)  (6) 

dx-»i 

dy->j 


This  is  analogous  to  the  continuum  definition  where  ijh  approaches  a 
point  of  zero  volume  on  a  macroscopic  scale  but  is  composed  of  discrete 
atoms  at  a  microscopic  level. 

For  simplicity,  a  microscopic  stiffness  coefficient  k  is 
introduced  to  represent  the  average  resistance  of  particles  to  movement 
in  the  z  direction  such  that  the  vertical  component  of  force  on  a 
particle  in  ijh  is  kw .  The  magnitude  of  k  depends  on  the  packing,  the 
roughness  of  the  particles,  and  the  confining  pressure.  The  vertical 
force  in  ijh  is  F  =  Nkw,  where  N  is  the  number  of  particles  in  ijh.  The 
vertical  stress  in  ijh  is 


Nkw  _  Nkh  ez 
Cz  (ij)  (ij) 


(7) 


Since  ijh  is  the  smallest  possible  volume  of  soil  that  can  be  used  for 
the  definition  of  stress  and  strain,  it  serves  as  a  control  volume 
analogous  to  a  point  in  a  continuum  and  particles  enter  and  leave  ijh  as 
deformation  takes  place.  This  means  that  N  varies  with  deformation,  and 
Equation  (7)  is  nonlinear.  This  can  easily  be  seen  as  N  =  i jh/[ Vp(l+e)], 
where  Vp  is  the  average  volume  of  a  particle  and  e  is  the  void  ratio; 
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and  by  definition,  the  volumetric  strain  is  related  to  the  void  ratio  as 
£v  =  (eo  -  e)/(l+eo).  where  e0  i s  the  initial  void  ratio.  Therefore 


_ Lit! _ 

CVp  (1  +  e0)  (1  -  cv) ] 


(8) 


Consider  two  elements  of  soil  adjacent  to  each  other  but  separated 
by  a  surface  dydz  in  the  yz  plane  as  shown  in  Figure  1.  Element  1  with 
center  at  location  x  has  particles  with  average  change  in  vertical 
displacement  w,  and  the  element  with  center  at  x+dx  has  average  change 
in  vertical  displacement  w+(3to/3x)dx .  The  vertical  force  at  the  right 
face  of  element  1  is  Fi  =  Nikw.  and  the  vertical  force  at  the  left  face 
of  element  2  is  F2  =  N2  k  [w+(3w/3x)dx]  where  Ni  is  the  number  of 
particles  at  the  right  face  of  element  1  and  N2  is  the  number  of 

particles  at  the  left  face  of  element  2.  The  vertical  shear  stress  at 
the  interface  between  the  elements  is  (Fi  -  F2 )/ ( dydz ) .  or 


txz  = 


[  (Ni  -  N2)  kw  -  N2ik|j] 

— 


(9) 


From  Equation  (7).  the  derivative  is 


3w 

3x 


It  ( iik CTz ) 


(10) 


The  substitution  of  Equations  (7)  and  (10)  into  Equation  (9)  gives 

3az 

Txz  =  baz  -  HD 


where 


b 


Ni  ~  n2  i  .  w  ik  3  ,ij_. 
N  h  2jh3x  Nk 


(12) 


and 
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Figure  1.  Adjacent  Soil  Elements. 


(13) 


It  can  be  noted  at  once  that  in  the  case  of  a  normal  load  changes 

of  i.  j,  and  N  with  x  can  be  assumed  small;  hence  Ni  =  N2  and  Equation 

(12)  gives  b  =  0.  However,  if  there  is  an  x  component  of  the  load, 

lateral  changes  of  N  with  x  may  not  be  neglected. 

1.  Two  Dimensions 

The  substitution  of  Equation  (11)  into  Equation  (2)  gives 

the  diffusion  equation 


daz 

dz 


a  a 

_ax(Dax)  ax^0^ 


(14) 


Equation  (14)  is  identical  to  that  of  References  9  and  11, 
however,  here  the  diffusion  coefficient  is  derived  from  the  relationship 
of  stress  to  particle  displacement.  Equation  (14)  is  nonlinear  as  D  is 
a  variable  that  depends  on  the  size  ijh  of  the  representati ve  element  of 
soil;  and  this,  in  turn,  depends  on  the  magnitude  of  the  particle 
displacements  or  stress  intensity.  For  small  particle  di spl acements ,  a 
small  volume  of  soil  ijh  may  be  used  to  adequately  define  stress  and 
strain  while  for  large  particle  movements,  a  larger  volume  of  soil  is 
needed  to  make  these  definitions.  If  the  horizontal  displacements  are 
small  compared  to  the  vertical,  then  ij  can  be  treated  as  a  constant  and 
h  is  seen  to  be  decreasing  away  from  the  loaded  area.  In  this  case,  D 
increases  with  z  and  the  absolute  value  of  x.  This  is  apparently  the 
same  type  of  diffusion  that  occurs  in  elastic  material,  for  it  can 
easily  be  shown  by  back  substitution  that  with  D  =  (x2  +  z2)/(2z)  and 
b  =»  q i D/ ( 2  +  qix)  where  qi  is  the  ratio  of  horizontal  to  vertical 

component  of  the  load.  Equation  (14)  gives  the  elastic  solution  for  a 
line  load  at  the  origin  of  coordinates. 

If  0  and  b  are  known,  then  with  known  boundary  conditions 
Equation  (14)  can  be  solved  for  the  vertical  stresses.  The  substitution 
of  Equation  (11)  into  Equation  (1)  gives  the  normal  horizontal  stress  as 
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0x  =  (D<T2)  + 


r°°  3 


—  (bGZ)dx 


(15) 


The  boundary  conditions  crx  =  az  =  0  at  x  =  -«>  gives 

0  f  OQ 

^  bazdx  =  0.  which  is  satisfied  only  if  b  is  a  constant.  Also 

from  Equation  (11)  at  x  =  0  and  z  =  0  it  is  known  that  3oz/3x  =  0. 

=  the  vertical  load,  and  ixz  =  the  horizontal  load;  hence  b  =  Qi. 
Therefore,  Equations  (11)  and  (15)  become 


3az 

txz  =  Ql<5z  ~  D"^~ 


(16) 


and 


<*X  =  0z’^crz)  +  dl  Txz 


(17) 


2.  Three  Dimensions 

In  three  dimensions,  the  equilibrium  equations  are 


3ax 

<kyx 

3tzx 

(18) 

3x 

+  3y  + 

~l 

II 

o 

3xXy 

3x 

3av 

+  if + 

^=0 

(19) 

dtxz 

dx 

+  %LL  + 
3y 

(20) 

txy 

=  Tyx.  tyz 

=  izy.  "txz  =izx 

(21) 

Development  similar  to  that  in  the  derivation  of  Equation  (11)  says  that 
changes  in  vertical  particle  displacement  across  the  xz  plane  produce 
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daz 

*xz  =  t>xz<*z  '  Dxz"^~ 


(22) 


and  changes  in  vertical  particle  displacement  across  the  yz  plane 
produce 


tyz  ~  tlyzCTz 


D 


daz 

y2“ 


(23) 


In  like  manner,  the  shear  stress  ixy  is  created  by  changes  in  particle 
displacement  in  the  y  direction  across  the  xy  plane  which  leads  to 

az 

*xy  ~  Dxy^z  '  Dxy^"  (24) 


The  subscripts  on  b  and  0  refer  to  b  and  D  in  the  respective  planes. 

If  the  inclined  load  has  X  and  Z  components  only,  then 
bxz  =  q i .  byz  =  0.  and  bXy  =  0.  Substituting  this  into  Equation  (20) 

gives  the  three-dimensional  diffusion  equation 


da  z 
Jz~ 


(25) 


Hence,  in  theory,  Equations  (18)  through  (25)  can  be  solved  for  all  of 
the  stresses. 


0.  SOLUTION  FOR  STRESSES  UNDER  INCLINED  LOADS 

1.  Two  Dimensional 

The  nature  of  D  in  Equation  (13)  is  unknown;  however,  as 
stated  above,  for  the  two  dimensional  case  D  is  a  function  of  i  and  h 
where  i  and  h  are  the  horizontal  and  vertical  components  of  the 
representati ve  element  of  soil.  The  magnitudes  of  i  and  h  depend  on  the 
magnitudes  of  the  particle  displacements  in  the  x  and  z  directions, 
respectively.  Directly  under  the  center  of  the  load  i  and  h  vary  with 
depth  only.  This  means  that  in  the  region  under  the  load  D  can  be 
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approximated  by  D(z),  a  function  of  z  only.  In  this  case,  one  can  make 
the  transformation 


W(z) 


D( z  )dz 


(26) 


and  the  solution  of  Equation  (14)  for  a  line  load  with  vertical 
component  P and  horizontal  component  qiP  is 


P 

°z  [4rtW(z)]i'2  exp 


( x- q j z  )2> 
4W(z )  , 


(27) 


The  shear  stress  can  then  be  found  from  Equation  (16)  as 


*zx  = 


(x-q1z)D(z)^ 

WT)  j  az 


(28) 


The  horizontal  stress  can  also  be  solved  from  Equation  (17)  as 


crx  = 


jo  *  (z) 


P2(  z  )  +  D2 (  z ) 
2W(z)  4W2(2) 


2qxW( z ) 


^2 1 


x  +  -mrl 


(29) 


where 


O’  (z) 


3D(z) 

3z 


(30) 


One  remarkable  property  noted  by  Golden  (1984)  is  that,  although 
nonrealistic  horizontal  tensile  stresses  are  predicted  by  elastic  theory 
if  horizontal  forces  are  present,  they  can  be  avoided  by  this  theory  if 
D ( z )  is  such  that  D ’ ( z )  2  D2 ( z ) / 2W ( z ) .  Another  observation  is  that  for 
a  purely  horizontal  force  the  elastic  theory  predicts  stresses  in 
granular  material,  while  in  these  equations  predict  zero  stresses. 
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2.  Three  Dimensions 

In  three  dimensions  with  Dxz  =  Dyz  =  D,  the  solution  of 

Equation  (25)  for  a  point  load  with  vertical  component  Qand  horizontal 
component  Q q i  in  the  x  direction  is 


oz 


Q 

4kW(z) 


exp 


i 


(x-q^)2  +  y2 
4W( z )  ' 


(31) 


It  can  also  be  shown  (Golden,  1986)  that  without  knowledge 
of  Dxy,  the  other  stresses  can  be  described  in  terms  of  az  by  the 

following  equations 


txz  =  q  iCFz 


ckr-i 


ax 


(32) 


tyz  =  'O' 


_3ctz 

w 


(33) 


txy 


3x3y  J 


Oz 


(34) 


ax  = 


Qi 


2  . 


+ D' + 


OZ 


(35) 


and 


ay 


+ 


where  D'  is  described  in  Equation  (30). 


(36) 
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SECTION  III 

STRAIN  ACCUMULATION  IN  GRANULAR  MATERIALS 
UNDER  STATIC  AND  REPEATED  LOADS 


A.  INTRODUCTION 

The  stress-strain  behavior  of  granular  soils  is  created  primarily 
b;  individual  particle  movements  to  form  a  denser  matrix  rather  than 
elastic  compression  of  the  particles.  As  a  result,  the  experimental 
curve  is  always  nonlinear  and  depends  on  the  loading  conditions.  In 
general,  the  shape  of  the  stress- strai n  curve  is  concave  towards  the 
strain  axis  in  triaxial  compression,  and  concave  towards  the  stress  axis 
in  one  dimensional  (uniaxial  strain)  and  isotropic  compression. 
Further,  it  is  noted  that  in  triaxial  conditions  some  sands  show  a 
distinct  peak  stress  while  others  do  not.  As  a  result,  it  is  difficult 
to  model  the  behavior  of  this  material  by  conventional  el asto- pi asti c 
theories,  and  recourse  is  usually  taken  to  empirical  methods  like  the 
hyperbolic  model  for  matching  triaxial  test  results  with  no  distinct 
peak  stress  (Oesai  and  Siriwardane,  1984)  or  numerical  curve  fitting 
techniques  (Desai.  1971). 

In  this  section,  traditional  approaches  to  stress-strain  modeling 
in  granular  media  is  identified  and  a  new.  more  rational,  approach  is 
introduced.  This  new  approach  is  particulate  and  is  in  agreement  with 
the  stress  transmission  theory  of  the  previous  section.  That  is.  the 
deformation  is  considered  to  be  the  result  of  particle  movement  rather 
than  deformation  of  particles.  The  result  is  a  general  model  that 
derives  the  stress - stra i n  response  for  one  dimensional,  isotropic,  and 
triaxial  loading  conditions  as  special  cases. 

B.  TRADITIONAL  APPROACH  TO  STRESS - STRAI N  MODELING  OF  GRANULAR 

MATERIAL 

The  experimentally-observed  stress - strai n  curves  of  granular  soils 
are  always  nonlinear.  The  shape  of  the  curve  also  differs  based  on  the 
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boundary  conditions  imposed  in  the  radial  directions.  As  a  result, 
traditional  elastic  and  elasto-plastic  theories  have  failed  to 
adequately  model  the  behavior  of  this  material.  So  far,  no  general 
model  exists  that  accounts  for  all  imposed  boundary  conditions  and 
traditional  methods  consist  of  applying  curve- fi tti ng  techniques  to 
curves  obtained  under  specific  conditions. 

Of  primary  importance  in  settlement  calculations  is  the  need  to 
model  the  triaxial  compression  condition.  One  simple  empirical  method 
proposed  by  Konder  and  Zelasko  (1963)  is  the  widely-used  hyperbolic 
moael  .  This  model  appears  to  be  a  natural  fit  to  soils  that  do  not 
exhibit  a  distinct  peak  stress.  For  soils  that  exhibit  a  distinct  peak 
stress,  a  modification  of  the  Ramberg-Osgood  empirical  model  used  for 
dynamic  loading  is  proposed  for  static  conditions  (Desai  and 
Siriwardane,  1984).  This  is  based  on  fitting  a  curve  to  the  initial 
tangent  modulus,  the  modulus  cf  the  plastic  zone,  the  yield  stress,  and 
a  parameter  defining  the  order  of  the  curve.  For  order  one,  this  curve 
becomes  a  hyperbola.  One  other  approach  suggested  is  to  treat  the 
tangent  shear  modulus  and  tangent  bulk  modulus  as  variables.  These  are 
usually  taken  as  linear  functions  of  the  octahedral  normal  and  shear 
stresses  (Nelson  and  Barron,  1971). 

The  difficulty  in  modeling  the  nonlinear  stress-strain  behavior  of 
sand  under  load  has  led  many  investigators  to  propose  numerical  curve¬ 
fitting  techniques.  One  of  the  most  popular  of  these  is  the  piecewise 
linear  method.  Here,  the  nonlinear  experimental  curve  is  divided  into 
pieces  of  linear  elastic  sections  for  numerical  analysis.  Very  often 
these  are  the  incremental  Hooke’s  law  or  the  hypoelastic  law  (Desai  and 
Wu,  1976).  Another  numerical  method  is  the  use  of  spline  functions  to 
fit  experimental ly-observed  curves.  These  are  functions  that  use  the 
data  to  provide  an  analytic  curve  similar  to  the  graphical  process  of 
using  a  French  curve  (Desai,  1971).  These  also  require  intense 
numerical  procedures,  and  the  data  must  be  presented  in  a  smooth  form 
and  not  scattered  as  observed  experimentally. 
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c. 


NEW  APPROACH  TO  STRESS - STRAl N  MODELING  OF  GRANULAR  MATERIAL 


1.  General  Theory 

Granular  media  is  composed  of  voids  and  particles. 
Deformation  is  caused  primarily  by  particle  movements.  As  the  particle 
movements  are  discrete,  their  derivative  at  a  point  does  not  exist,  and 
the  strain  in  this  media  cannot  be  defined  at  a  point.  Strain  can  only 
be  defined  with  respect  to  an  elemental  volume  of  soil  with  enough 
particles  so  that  the  relative  movements  of  the  particles  in  the 
elemental  volume  can  produce  deformation.  Like  the  strain,  the  stress 
in  granular  soils  does  not  exist  in  a  void  and  should  not  be  described 
at  a  point.  The  stress  correspondi ng  to  the  strain  in  the  elemental 
volume  is  the  result  of  forces  on  the  particles  in  the  elemental  volume 
and,  like  the  strain,  can  only  be  defined  with  respect  to  the  elemental 
volume.  These  definitions  of  strain  and  stress  are  given  analytically 
as  Equations  (5)  and  (6)  of  Section  II. 

The  substitution  of  Equation  (5)  into  Equation  (6)  gave  the 
stress -strai n  relationship  of  Equation  (7).  Referring  to  the  z 
direction  as  the  axial  direction  and  letting  subscript  "a"  refer  to  the 
axial  direction,  the  stress-strain  relationship  in  the  axial  direction 
as  obtained  from  Equation  (7)  is 


Nkh  ea 
CTa  =  (i  j) 


(37) 


where  N  is  a  function  of  the  volumetric  strain  as  shown  in  Equation  (8). 
The  substitution  of  Equation  (8)  into  Equation  (37)  gives  the  stress- 
strain  relationship  in  the  axial  direction  as 


_  £q  £a 

°a  "  (l-jJev) 


(38) 


where 
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(39) 


=  h2  k 

°  [vp  ( l+e0)  ] 


Here,  (3  =  l/evi  and  evi  is  the  maximum  value  obtainable  by  ev .  This 
latter  term  was  added  since  ev  =  1  i s  unattainable  under  conventional 
loads . 

The  derivative  of  Equation  (38)  gives  the  slope  of  the 
stress-strain  curve  as 


(40) 


It  is  apparent  from  this  that  at  ea  =  ev  =  0,  the  slope  is  E0 . 
Therefore.  E0  is  the  initial  tangent  modulus  of  the  soil,  and  knowledge 
of  h  and  k  is  not  necessary  if  E0  can  be  measured. 

2.  Application  To  One  Dimensional  And  Isotropic  Loading 
Conditions 


In  one  dimensional  compression  of  soils  no  lateral  strains 
are  allowed,  and  the  volumetric  strain  is  equal  to  the  axial  strain. 
That  is  ev  =  ea  and  evi  =  eai_.  where  eaL  is  the  asymptotic  axial  strain. 

Therefore,  in  this  case,  Equation  (38)  becomes 


Eq  ea 


L 


(41) 


In  isotropic  loading  conditions,  the  strains  are  equal  in  all 
directions.  Hence.  ev  =  3ea  and  evi  =  3eaL,  and  Equation  (41)  also 

applies  for  this  case. 
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It  should  be  noted  that  Equation  (41)  can  be  written  as 


cra  =  £c 


(£a  +  —  + 
£aL 


.  .  ) 


(42) 


and 


2~— •  [exp(2  ea/eaL  -  1)  -  1  ]  =  E0  (ea  +  ea/EaL  +  ...)  (43) 

Therefore  letting  04  =  2/eat_  and  C  =  E0/a!  gives  the  alternate 
expression 

aa  =  C  expict^)  -  C  (44) 

Equation  (44)  is  the  same  as  that  derived  by  the  hypoelastic  analysis 
for  these  loading  conditions  (Desai  and  Siriwardane.  1984). 

3.  Application  To  Triaxial  Loading  Conditions 

In  elastic  material  under  triaxial  compression  the  radial 
strain,  £r,  is  proportional  to  the  axial  strain,  and  the  proportional 
constant.  Poisson's  ratio,  is  determined  empirically.  In  effect,  if  vt 
=  -der/dea  and  vs  =  *er/ea-  then  in  elastic  material  vt  =  vs.  However, 
in  triaxial  compression  of  sands,  the  relationship  of  radial  strain  to 
axial  strain  is  nonlinear.  To  represent  this  nonlinear  condition  one 
can  let  vs  -  Vt  =  D2.  where  Dzis  an  empirical  constant  representing  the 
average  difference  between  v$  and  vt.  If  D2  =  0  the  relationship  is 
linear,  and  the  nonlinearity  increases  as  the  magnitude  of  O2  increases. 
The  definition  vs  -  -er/ea  gives 

dvs  2  Dt 

3^-  (er  t  Ca  Vt/Ea)  -  -  _  (45) 
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Also  since  ev  =  ea+2er  then  Vs  =  (£a  '  ev)/(2ea).  which  shows  that 


dv$  d ( £y /ea ) 

aJT  "  "  °’5  — rer~ 


(46) 


Equating  Equations  (45)  and  (46)  gives  d(£v/£a)/d£a  =  2D2/ea.  which  has 
sol uti on 

ev  =  8263  +  2D2£alnea  (47) 

where  B2  is  an  integration  constant.  It  should  be  noted  that  xlnx 
approaches  zero  as  x  approaches  zero.  For  example.  (0 .OODlnO.OOl  =  - 
0.007  and  ( 0. 0001 ) InO . 0001  =  -0.0009. 

The  substitution  of  Equation  (47)  into  Equation  (38)  yields 

_ £q  £a _ 

°a  (1  -  pB2ea  -  2p02ealnea)  C 

Letting  cra  =  au  and  ea  =  eu  at  maximum  stress.  Equation  (48)  gives 
PB2  =  l/eu  -  E0/au  *  2p02lneu.  Also  setting  d<ra/dea  *  0  at  maximum  stress 
given  2PD2  =  -l/eu.  The  substitution  of  these  into  Equation  (48)  yields 
the  general  relationship 

aa  =  T - - - T  (49) 

[_a  +  a(ea/eu)lnea  +  bea  J 

where  a  =  1/E0  and  b  =  l/ou  -  ( l+lneu)/(  E0eu)  •  In  soils  with  no  distinct 
peak  stress  such  as  loose  sands  and  sands  under  high  pressure.  eu 
approaches  infinity,  therefore.  Equation  (49)  reduces  to  the  hyperbolic 
model 


°a  (  a  +  bea  ) 
where  a  =  1/E0  and  b  =  1/ctu. 


(50) 
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The  constants  82  and  O2  in  Equation  (47)  can  also  be  evaluated  in 
terms  of  the  values  at  maximum  volumetric  strain.  If  ea  =  e0  at 
dev/dea  =  0.  then  82  =  -  2D2(  l  +  lne0 ) .  Also,  since  at  ea  =  e0  the 
volumetric  strain  ev  =  evm.  the  maximum  volumetric  strain.  Equation  (47) 

becomes 

£y  =  ea(£vm/eo  1  Cl  In  (e0/ea) ]  (51) 

This  equation  shows  that  ev  =  0  at  the  two  points  ea  =  0  and  ea  =  ei. 
The  value  of  ey  is  obtained  from  Equation  (51)  as  ei  =  2 . 7 18e0 .  Further, 
comparison  of  these  values  of  B2  and  D2with  those  found  in 
Equation  (48)  gives  (5  =  £0/(euevm)  and  eu  satisfying  E0eu  -  auln(eu/e0). 

0.  PERMANENT  AXIAL  STRAIN  ACCUMULATION  IN  GRANULAR  MATERIALS 

1.  Static  Loading 

Equation  (49)  gives  the  axial  stress - strai n  relationship 
under  triaxial  conditions.  In  this  equation,  ea  is  the  total  strain 
accumulated  under  the  static  stress  increment  ct3.  Upon  removal  of  the 
load,  the  soil  rebounds.  The  amount  of  rebound  is  the  elastic  portion 
of  the  total  strain  and  can  be  written  as 


where  Er  is  the  resilient  modulus  of  the  soil.  The  permanent  strain 
accumulated  due  to  one  application  of  oa  is 

ea  =  ea  -  ea  =  ^  h(oa)  -  ^  <Ja  (53) 

where  h(oa)  is  obtained  from  Equation  (49)  by  solving  for  ea  in  terms  of 
oa .  For  the  hyperbolic  case  it  is  obtained  from  Equation  (50)  as 
h(oa)  =  aaa/(l-boa).  Equation  (33)  may  be  written  in  the  more  compact 
form: 
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(54) 


where 

f(aa)  =  h(aa)  '  ki<ja  (55) 

and  kj  =  E0/Er. 

2.  Repeated  Loading 

At  every  load  application  there  is  some  readjustment  of  the 
grains  in  a  soil  sample.  This  rearrangement  becomes  less  pronounced  as 
the  number  of  load  applications  increases  because  the  soil  becomes  more 
packed.  Therefore,  the  rate  of  change  of  initial  stiffness  at  any  load 
application  depends  on  the  particle  arrangement  (or  on  the  stiffness)  at 
that  load  application  and  decreases  with  increasing  load  applications. 
The  simplest  way  to  represent  this  behavior  is  by  a  power  function 


KiE 


01 


(56) 


where  Eoi  is  the  initial  target  modulus  at  the  ith  cycle  of  loading.  Ki 

is  a  proportional  constant,  and  n  is  a  parameter  reflecting  the 
dependence  of  the  rate  of  change  of  Eoi  on  Eoi.  The  solution  of 
Equation  (56)  gives 


Eoid^oi 


En+1 

^01 

n+1 


Ai 


(57) 


where  A]  is  an  integration  constant.  Solving  for  E0i  gives 

r  . 

Eoi  =  CAi  +  B] 


(58) 


where 
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A  =  ( n+1 ) Ki ,  B  =  - Aj(  n+1 ) ,  m  =  ~ 


(59) 


The  permanent  strain  due  cycle  i  is  then  obtained  from 
Equation  (54)  as 


e 


P 

ai 


=  rr  fi(aa) 
t-oi 


(60) 


where  fi(aa)  is  the  relation  defined  by  Equation  (55)  for  the  ith  load 
appl i cati on . 

The  permanent  strain  accumulated  in  N  cycles  is  then 


p 

eaN  ~ 


f i ( aa )  [Ai  +  B]'m  di 


J  o 


(61) 


For  soils  where  the  hyperbolic  model  presents  a  good  fit.  it  can  be 
shown  that  f-j(<ra)  =  f(cra)  if  <ra  is  constant,  that  is,  it  is  independent 

of  load  cycle.  This  is  because  for  these  soils,  the  assymptotic  stress 
ctu  and  the  ratio  ki  =  E0/Er  are  fairly  constant  during  load 
applications.  For  this  case.  f-j(cra)  can  be  taken  out  of  the  integral, 
and  the  result  of  the  integration  of  Equation  (61)  depends  on  whether  m 
i s  equal  to  one  or  not . 

a .  Case  Where  m  =  1 

In  this  case,  n  =  0  or  the  change  of  Eoi  where  load 
application  is  constant.  In  this  case,  Equation  (61)  becomes 


P 

eaN 


f  (aa) 


CN  di 

f  (aa) 

m( 

>  +  £) 

■  in£] 

Ai+B 

A 

J  0 

(62) 


From  Equation  (59)  for  n  =  0,  it  is  seen  that  B/A  =  -Ai/K],  and  from 
Equation  (57)  for  i  =  1  and  i  =  2,  one  finds  Ki  =  E02  -  E01  and 
At  =  E02  •  2E0i.  This  means  that 
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(63) 


B_  =  2Eqi  -  Eo2 
A  E02  '  E0i 


In  general,  E02  is  greater  than  E01 :  hence  B/A  is  small  compared  to  N. 
This  means  that  as  an  approximation,  the  permanent  axial  strain 
accumulated  is  obtained  from  Equation  (62)  as 

p 

eaN  =  ai  +  bi  In  N  (64) 


Where  bi  =  f(oa)/A  and  ai  =  ( 1 / A ) f ( <ya )ln(B/A)  are  constants.  It  is  also 

□ 

seen  from  Equation  (64)  that  ai  =  eai  as  defined  in  Equation  (60). 
Equation  (64)  is  that  proposed  by  Lentz  and  Baladi  (1981)  based  on  its 
goodness  of  fit  to  their  experimental  data, 
b.  Case  Where  m  t  1 

In  this  case,  Equation  (61)  becomes 


eaN  =  f(<*a) 


[Ai  +  B] 'm  di 


A1'1”  f(ca) 
1-m 


1-m 

1-mi 

[n+a] 

'  1 

(a) 

I 

(65) 


Again  from  Equation  (59)  it  is  seen  that  B/A  =  - Aj / K 1 ,  and  from 
Equation  (57),  that  (n+l)Kl  -  E^1-  EflJ^and  (n+l)Ai  =  E^1  -  2  EqI  1 . 

From  this.  B/A  is  obtained  as 

2C11-  EnoT 

7-  .rH-l  "vr  (66) 

Eol  '  eo2 


Since  ^02  is  greater  than  E01,  it  is  seen  that  B/A  i$ 
small  compared  to  N,  and  Equation  (65)  can  be  approximated  by 


P 

EaN 


=  Bi  N 


-m 


(67) 
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where  Bj  =  a1  - nrif (<ja)/(  1-m) .  This  equation  is  the  same  as  that  suggested 
by  Khedr  (1986)  and  Diyaljee  and  Raymond  (1982)  based  on  their 
experimental  work.  It  should  also  be  noted  that  from  Equation  (67) 

p 

Bj  =  eai  as  defined  in  Equation  (60). 

E.  PERMANENT  RADIAL  STRAIN  ACCUMULATION  IN  GRANULAR  MATERIALS 


The  total  radial  strain  during  static  loading  is  er  =  (e„  -  ea)/2 
where  ev  is  given  by  Equation  (51).  The  elastic  radial  strain  is 

e  e  v  ,co. 

£ r  =  - vea  =  *7-  ca  (68) 

Lr 


where  vis  the  elastic  Poisson  ratio  during  unloading.  This  means  that 
the  permanent  radial  strain  is 

e?  =  er  -£p  (69) 


The  cumulative  radial  strain  for  N  cycles  of  loading  is  then 


||^eai  [l-  ln(e0/eai)] 


di  (70) 


where 


£ai  ~  r  .  h(aa )  (71) 

Loi 

and  E01-  is  given  by  Equation  (58).  The  integration  presented  by 
Equation  (70)  is  clumsy.  An  approximate  more  straightforward  estimate 
can  be  evaluated  using  the  formula  presented  by  Chang  and  Whitman 
(1988).  In  their  evaluation,  they  found  that 
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(72) 


=  M2  -n2 

P  2t| 

Y 

where  +  2  is  the  permanent  volumetric  strain,  =  2/3(e^-6p) 

is  the  permanent  shear  strain.  Also,  M  is  the  ratio  of  mean  and 
deviatoric  stress  at  ea  =  £o .  that  is  at  minimum  volumetric  strain,  and 
r|  is  the  ratio  of  mean  and  deviatoric  stress  at  ea  =  0.  or  at  the  start 
of  loading.  The  mean  stress  is  (oa  +  2 crr)/3,  and  the  deviatoric  stress 
is  oa  '  <*r.  From  Equation  (72).  the  permanent  radial  strain  can  be 
solved  in  terms  of  the  permanent  axial  strain  as 

A  ,  M?  -!L:  33.  »  (73) 

M2  -  t)2  +  6t\ 


The  permanent  radial  strain  accumulated  after  N  cycles  is 


p  _  m2  -  n  -  3n  p 

er"  '  M2  •  *  +  6n  £a" 


(74) 


where £gN  is  the  permanent  axial  strain  accumulated  after  N  cycles  as 
given  in  Equation  (64)  or  (67). 
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SECTION  IV 

VALIDATION  OF  PART  I ULATE  THEORY  FOR  GRANULAR  MATERIAL 


A.  INTRODUCTION 

In  this  chapter,  the  assumptions  made  in  the  theory  leading  up  to 
the  stress  continuity  equation  of  particulate  media  (Equation  11).  the 
constitutive  equation  of  particulate  media  (Equation  38),  and  the  strain 
growth  equation  of  particulate  media  (Equation  61)  is  verified.  The 
method  of  validation  is  by  the  comparison  of  these  analytical 
expressions  with  experimental  observations  found  in  the  literature. 

To  validate  the  stress  continuity  concept  for  granular  media,  the 
stresses  predicted  by  the  particulate  theory  developed  in  Section  II  is 
evaluated  against  measured  stresses  in  particulate  media.  The  results 
are  further  compared  with  the  predictions  of  the  linear  elastic  theory. 
The  stress-strain  theory  for  granular  media  developed  in  Section  III  is 
also  compared  with  measured  stress-strain  data  on  granular  media. 
First,  the  stress-strain  responses  under  static  isotropic  and  static 
triaxial  compression  are  examined.  Next,  the  theory  on  permanent  strain 
accumulation  in  granular  material  under  repeated  loading  is  compared 
with  experimental  observations  using  different  load  magnitudes. 

B.  VALIDATION  OF  PARTICULATE  THEORY  OF  STRESS  TRANSFER 

It  was  already  shown  in  Section  II  that  the  particulate  theory  of 
stress  transfer  in  granular  media  does  not  predict  horizontal  tensile 
stresses  for  inclined  loads  and  predicts  no  stress  transfer  for  purely 
horizontal  loads.  These  two  phenomena  are  unique  to  granular  material 
and  cannot  be  avoided  by  any  theory  that  imposes  the  requirement  for  the 
existence  of  the  second  derivatives  of  strains  in  the  media.  In  this 
section,  the  comparison  of  the  theory  for  the  prediction  of  vertical 
stresses  under  a  loaded  area  with  observed  experimental  results  will  be 
made . 
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Morgan  and  Gerrard  (1981)  presented  experimental  results  of 
several  investigators  on  the  stress  distribution  in  sands  under  load. 
Fairly  compl ete  i nformat i on  was  obtainable  from  their  report  on  the 
results  of  the  Waterways  Experiment  Station  (WES)  tests  and  the  tests  at 
Melbourne  University  (Ml  and  M2).  These  tests  were  performed  using 
uniform  vertical  loads  (q)  over  circular  areas  of  radius  "a"  and  the 
results  included  the  distribution  of  vertical  stress  with  depth  directly 
under  the  center  of  the  areas  as  shown  in  Figure  2. 

Directly  under  the  center  of  a  circular  area  of  radius  "a",  the 
elastic  solution  predicts  from  a  vertical  uniform  load  q 

f-  1.[(./Z)»*  l]'3'2  (75) 

This  is  shown  as  the  dashed  line  of  Figure  2. 

For  a  vertical  uniform  distributed  load  of  intensity  q  over  a 
circular  area  of  radius  "a",  with  coordinate  origin  taken  at  the  center 
of  the  loaded  area.  Equation  (31)  gives 

°z  =  4kH(2")  fZ"  PexP  [  ‘  4wT77  (r2  +  p2  '2rpcos0)  ]  rdrdG  ( 76 ) 

J  o  J  o 

where  p2  =  x2  +  y2.  That  is.  in  Equation  (31)  let  qi  =  0,  replace  Q  by 
qrdrd0,  and  employ  the  law  of  cosines  (Figure  3).  Directly  under  the 
center  of  the  area,  p  =  0.  and  the  integration  of  Equation  (76)  gives 

Y  =  1-exp  {  -a2/[4W(z)]  }  (77) 

Reference  11  showed  that  the  assumption  that  D(z)  is  a  linear 
function  of  z  satisfies  all  granular  soil  requirements.  From 
Equation  (26).  this  means  that  W(z)  can  be  expressed  as 

Wi|l  =  dt  +  d2(  |)  (78) 
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A. 
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The  substitution  of  Equation  (78)  into  Equation  (77)  gives  the 
transformed  equation 


-  (a/z)2 
41 n( l-oz/q ) 


=  di  +  d2  (  |-  ) 


(79) 


This  equation  is  suitable  for  linear  regression  on  the  data  of  Figure  2. 
This  linear  regression  produced  the  solid  lines  in  Figure  2.  For  each 
soil  the  fit  is  remarkable,  indicating  that  the  form  of  Equation  (78)  is 
appropriate  for  this  type  of  loading.  The  regression  revealed  that 
di  =  0.078,  d2  =  0.217  for  the  Ml  data,  di  =  0.032,  d2  =  0.157  for  the 
M2  data,  and  di  =  0.101,  d2  =  0.057  for  the  WES  data. 

Figure  2  can  also  give  an  indication  of  the  relative  lateral 
spread  of  the  distribution.  To  do  this,  the  approximate  method  of 
analysis  for  vertical  stress  directly  below  the  loaded  area  is 
employed.  This  says  for  p  =  0 

oz  =  - - -  (80) 

n  (a  +  ztan  a)2 

Where  tan  a  is  a  measure  of  the  magnitude  lateral  spread  with  depth.  A 
value  of  a  =  30  degrees  estimated  from  the  theory  of  elasticity,  is 
usually  assumed  (Dunn,  Anderson,  and  Kiefer,  1980).  From  Equation  (80), 
it  is  found  that 


(81) 


The  best  fit  value  for  the  angle  a  from  the  data  of  Figure  2  is 
a  =  26.01  degrees  for  MI,  a  =  19.33  degrees  for  Mil,  a  =  23.81  degrees 
for  WES  and  a  =  29.31  degrees  for  the  elastic  prediction.  The  larger 
angle  of  a  for  the  elastic  prediction  indicates  that  the  elastic 
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solution  predicts  a  wider  lateral  spread  of  the  vertical  pressure 
distribution.  The  method  of  compaction  is  seen  to  be  a  determining 
factor  in  the  spread,  as  the  MI  sand  was  identical  to  the  Mil  sand 
except  that  the  former  was  compacted  by  a  vibrating  plate  and  the  latter 
by  pluvial  compaction. 

C.  VALIDATION  OF  PARTICULATE  THEORY  OF  STRESS - STRAI N  RESPONSE 

It  was  shown  in  Section  III  that  Equation  (38)  is  a  general 
expression  for  the  relationship  of  axial  stress  to  axial  strain.  For 
the  special  conditions  of  one  dimensional  and  isotropic  compression, 
this  equation  reduced  to  Equation  (41)  and,  for  the  case  of  isotropic 
compression,  to  Equation  (49). 

The  validation  of  Equation  (41)  was  made  using  data  for  isotropic 
loading  conditions  on  two  samples  of  McCormic  Ranch  sand  presented  by 
Desai  and  Siriwardane  (1984,  pg.  193)  and  shown  in  Figure  4.  In  both 
cases,  the  regression  of  Equation  (41)  on  the  data  revealed  99.9  percent 
correlation.  The  regression  curve  is  shown  as  the  dashed  lines  in 
Figure  4. 

The  validation  of  Equation  (49)  was  also  made  on  data  for  a  medium 
dense  sand  presented  by  Desai  and  Siriwardane  (1984,  pg  177).  The 
measured  stress-strain  response  of  this  sand  at  two  different  confining 
pressures  (03)  are  shown  in  Figure  5.  Again,  the  fit  of  the  equation 

was  excellent  and  a  regression  of  Equation  (49)  on  the  data  revealed 
correlation  coefficients  of  99.8  percent.  The  predicted  curves  are 
shown  as  the  dashed  line  in  Figure  5. 

Equation  (49)  is  based  on  the  nonlinear  relationship  of  the 
volumetric  strain  to  the  axial  strain  derived  in  Equation  (47).  This 
equation  was  also  evaluated  using  the  volumetric  versus  axial  strain 
measurements  of  the  sand  in  Figure  5  taken  at  the  two  confining 
pressures.  This  data  is  shown  in  Figure  6.  The  regression  of 
Equation  (47)  on  this  data  is  shown  as  the  dashed  line  in  the  figure. 
The  correlation  coefficient  was  99  percent  for  each  case. 
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Figure  4.  Stress-Strain  Behavior  For  Isotropic  Loading 

(Data  From  Reference  6,  pl93,  1  psi  =  6.895  kPa). 
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Equation  (49)  with  £  =  0.021 

a  =  0.0012/  a, 
b  =  0.3601/ 


r/ 


Equation  (49)  with  £ 
a  =  0.0035/  (j-\ 
b  =  0.2418/  <r3 


X  Data  Points  For  =  0.49  kg/cm, 

A  Data  Points  For  <73  =  9.79  kg/cni 
—  Regression  Curves  J99.8%  Correlation 


Figure  5.  Stress-Strain  Behavior  For  Triaxial  Compression 
(Data  From  Reference  6,  p 177 ) . 


o  0.01  0.02  0.03  3«a 


Figure  6.  Relationship  Of  Volumetric  And  Axial  Strains 
(Data  From  Reference  6,  p 177 ) . 
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D. 


VALIDATION  OF  THE  PARTICULATE  THEORY  OF  STRAIN  ACCUMULATION  UNDER 
REPEATED  LOADING 


The  nature  of  strain  accumulation  in  granular  material  is 
presented  in  the  general  form  of  Equation  (61).  This  equation  is  a 
function  of  a  parameter  m  that  reflects  not  only  the  nature  of  the 
material  Put  also  the  magnituae  of  the  applied  load.  To  see  this, 
Equation  (39)  can  be  written  as 


h2k( l-n0) 


(82) 


wnere  n0  =  e0/(l+e0)  is  the  initial  porosity  of  the  soil.  For  small 
loads,  the  soil  compacts  in  fairly  equal  increments,  and  changes  in  the 
porosity  are  of  equal  increments.  From  Equation  (82),  this  means  that 
for  small  loads,  the  initial  tangent  modulus  increases  linearly  with  the 
number  of  load  applications.  In  this  case.  Equation  (58)  says  that  m  is 
close  to  unity  and  Equation  (61)  gives  the  semi-log  relationship  of 
Equation  (64).  For  large  loads,  the  soil  compacts  in  more  uneven 
increments  with  larger  increments  occurring  at  the  early  cycles  of 
loading.  This  reveals  the  m  in  Equation  (58)  should  be  less  than  unity, 
and  Equation  (61)  becomes  the  log-log  expression  of  Equation  (67). 

To  validate  Equations  (64)  and  (67).  data  on  permanent  axial 
strain  accumulation  with  number  of  cycles  of  loading  for  five  load 
magnitudes  on  a  Dolomite  Ballast  presented  by  Diyaljee  and  Raymond 
(1982)  were  used.  These  experimental  results  are  plotted  on  a  log-log 
plot  in  Figure  7.  In  each  case.  Equation  (67)  represented  a  perfect 
fit,  with  m  =  0.87  for  the  four  smaller  loads  and  m  =  0.80  for  the 
largest  load. 

The  logic  behind  Equation  (64)  can  also  be  tested  using  the  data 
of  Figure  7.  To  do  this,  the  data  are  also  plotted  on  a  semi-log  scale 
as  shown  in  Figure  8.  The  data  with  m  closest  to  unity  should  plot  more 
linear  than  those  with  values  of  m  furthest  from  unity.  This  is  seen  in 
Figure  8  where  the  data  with  m  =  0.87  also  plotted  a  straight  line  on 
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Permanent  Axial  Strain  (%) 


Figure  8.  Semi-Log  Plot  Of  Permanent  Axial  Strain  Versus 
Number  Of  Load  Cycles  (Data  From  Reference  7). 
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this  semi-log  scale.  However,  the  case  of  the  heaviest  load,  with 
m  =  0.80,  shows  a  distinct  non-linear  plot  on  the  semi-log  scale. 
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SECTION  V 

EVALUATION  OF  MATERIAL  PARAMETERS  IN  PARTICULATE  THEORY 


A.  INTRODUCTION 

In  the  development  of  the  stress  distribution  theory  of  Section  II 
ana  the  stress-strain  and  strain  accumulation  theories  of  Section  III, 
several  material  parameters  were  identified.  These  parameters  are 
material  specific  and.  like  Young's  modulus  and  Poisson's  ratio  of  the 
theory  of  elasticity,  are  best  evaluated  from  experimental  observations. 
However,  it  must  first  be  determined  if  these  parameters  are  also 
functions  of  geometry  and  load.  If  they  are.  their  dependence  on  these 
quantities  should  be  identified  such  that  experimental  procedures  can  be 
scaled  to  observe  only  the  material  characteristics  of  the  parameters. 

In  this  section,  the  material  parameters  arising  in  the  stress 
distribution  theory,  the  stress - strai n  theory,  and  the  strain 
accumulation  theory  are  examined  and  methods  of  determining  them  are 
discussed.  In  particular,  their  relationship  to  other  parameters  such 
as  depth  in  the  soil,  confining  pressure  and  load  intensity  are 
i denti f i ed . 

8.  MATERIAL  CONSTANTS  IN  PARTICULATE  THEORY  OF  STRESS  TRANSFER 

The  theory  of  stress  transfer  in  granular  material  presented  in 
Equation  (14)  depends  on  two  constants  D  and  b.  It  is  shown  in 
Section  II  that  b  =  di.  the  ratio  of  the  horizontal  component  to  the 
vertical  component  of  the  applied  load.  The  parameter  D  is  the 
diffusion  coefficient  of  the  soil.  In  Equation  (13)  N2  is  a  fraction  of 
N;  therefore  D  is  dependent  only  cn  the  size  of  the  elemental  volume  ijh 
needed  to  adequately  define  stress  and  strain  in  the  material.  In 
particular,  D  is  proportional  to  the  square  of  the  horizontal  component, 
i2,  and  inversely  proportional  to  the  vertical  component  h.  For  small 
particle  displ  a  cement .  a  small  vol  ume  ijk  of  soil  may  be  used  to 
adequately  define  stress  and  strain,  while  for  large  particle 
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ci  sp :  acemerus ,  3  larger  o  1  urr.e  cr  aci:  :  •;  needed  to  make  these 
aefinitions.  As  a  result,  the  magnitudes  cf  1.  j,  ana  h  at  a  given 
location  are  dependent  on  tne  stress  level  at  that  location.  It  is 
shown  in  Section  IV  that  for  a  vertical  lead  on  a  circular  foundation  of 
raGius  "a",  me  coefficient  D  can  oe  approximated  by  a  linear  function 
of  oeptn  z  as 


3  =  2d i z  -r  ado  (83) 

where  dj  a  no  do  are  material  parameters  to  be  determined  from 

experiments.  Tne  linear  relationship  of  D  to  z  is  also  observed  by 

Golden  (1986)  and  Hill  and  Harr  ',1971).  Equation  (83)  indicates  that 
ado  is  a  term  reflecting  the  boundary  condition  at  z  =  0,  and  di  is  a 

parameter  reflecting  the  lateral  spread  of  the  vertical  stress 

di stri Duti on . 

As  shown  in  Section  IV.  both  di  and  d2  are  dependent  on  the  method 

of  compaction  of  the  granular  soil.  However,  an  accurate  determination 
of  the  stresses  below  the  loaded  area  can  be  made  if  di  and  d2  are 

determined  for  that  soil.  Since  both  di  and  d2  are  material  parameters, 
this  can  be  determined  by  laboratory  tests  for  small-scale  loads.  The 
linear  nature  of  Equation  (76)  in  terms  of  loads  indicates  that  the 
values  of  dj  and  d20btained  in  the  laboratory  can  be  used  for  field 
predictions.  The  determination  of  dj  and  d2  in  the  laboratory  requires 
the  measurement  of  vertical  pressure  with  depth  directly  below  the 
center  of  the  loaded  area  as  illustrated  by  Morgan  and  Gerrard  (1981). 
The  values  of  di  and  d2  are  then  evaluated  by  regression  of 

Equation  (79)  on  tne  data. 

C.  MATERIAL  CONSTANTS  IN  PARTICULATE  THEORY  OF  STRESS  -  STRAIN  RESPONSE 

The  general  stress - strai n  response  of  particular  media  is  given  by 
Equation  (38).  For  one  dimensional  and  isotropic  loading,  this  reduces 
to  Equation  (41)  and  for  triaxial  loading  to  Equation  (49).  Since  the 
triaxial  condition  is  applicable  directly  to  the  rutting  problem,  the 
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mature  or  t.ne  cons  tar.  is  '  n  Eauaticn  (49)  is  examined  here. 
Equation  ■* 49 )  can  oe  rewritten  as 


<?a 


_ -o£a _ 

l  -*■  ( £3/ -£a/£u)C^rl  ( £a^£u )  '  1  ] 


(84) 


wnere  e5  =  au7E0  is  the  strain  predicted  at  maximum  stress  if  the 
material  was  elastic  and  eu  is  the  strain  observed  at  maximum  stress. 
This  equation  is  hyperbolic  if  eu  is  infinite  and  is  linear  if  ee  =  eu 
and  finite  as  in  the  elastic  case.  The  three  material  parameters  E0. 
£e .  and  eu  will  be  examinee  separately  in  this  section. 

The  initial  soil  modulus  is  developed  in  Equation  (39)  as  a 
function  of  the  particle  size,  initial  void  ratio,  the  height  h  of  the 
minimal  volume  of  granular  soil  necessary  to  determine  stress  and 
strain,  and  a  parameter  k  representing  the  average  resistance  of  a 
particle  to  movement  in  the  z  direction.  The  parameters  h  and  k  are  not 
easily  determined.  Fortunately,  the  initial  modulus  can  be  obtained 
from  the  regression  of  Equations  (49)  or  (84)  on  measured  stress-strain 
data.  As  the  magnitude  of  k  depends  on  the  roughness  of  the  particles, 
the  packing,  and  the  confining  pressure,  the  dependence  of  the  initial 
soil  modulus  on  the  confining  pressure  is  apparent  from  Equation  (39). 
This  dependence  as  observed  experimentally  is  usually  expressed  as 
(Seed,  et  al  ,  1986;  Richart.  Hall  and  Woods,  1970  ) 


E0 


Ci 


(85) 


where  Ci  is  a  constant  and  a3  is  the  confining  pressure.  Figure  9  shows 
that  the  relationship  holds  for  the  medium  dense  sand  presented  by  Desai 
and  Si riwardane  ( 1984 ) . 

Figure  9  shows  that  ee  and  eu  are  also  proportional  to  the  square 
root  of  the  confining  pressure.  Since  the  confining  pressure  insitu  is 
the  product  of  the  coefficient  of  earth  pressure  at  rest,  the  effective 
unit  weight  and  the  depth,  the  proportional  relationships  of  E0,  ctu>  and 
eu  to  the  square  root  of  the  confining  pressure  allow  extrapolation  and 
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;  nterpol ation  using  the  experimental  results  of  only  one  confining 


pressure . 


material  constants  in  particulate  theory  of  strain  accumulation 


'he  accumulation  of  permanent  stress  after  N  cycles  cf  load 
applications  is  given  by  Equation  (67),  which  for  small  values  of  n  can 
ce  approximated  by  Equation  (64).  in  both  equations,  the  constant  can 
ce  obtained  oy  regression  of  the  equations  on  experimental  data, 
unf crtuna tel y ,  tnese  constants  are  functions  of  the  magnitude  of  the 
applied  load.  In  Equation  (67)  the  constant  Bi  is  the  permanent  strain 
accumulated  after  the  first  application  of  the  load  and  1-m  is  the  slope 

of  tne  log-log  transformation  of  the  equation.  In  Equation  (64),  the 
constant  at  is  tne  permanent  strain  accumulated  after  the  first 
application  of  the  load,  and  bj  is  the  slope  of  the  semi-log 
relationship. 

The  permanent  strain  accumulated  after  the  first  application  of 
the  load  is  given  in  Equation  (53).  However,  in  the  deviation  of 
Equation  (67).  it  was  seen  that  Bi  =  A1  'mf  ( aa )  /  ( 1  -m)  where  from 
Equation  (59)  m  =  l/(n+l)  and  A  =  (n+l)Ki,  and  from  Equation  (57)  it  is 
seen  that  (n+l)Ki  =  Eg^1-  E1^.1  These  developments  show  that  m  satisfies 

the  relationship 


1  /rru  1  /m 


“01 


"V 


m 


(86) 


The  value  of  aj  in  Equation  (64)  is  given  by  Equation  (53).  Also  in  the 
deviation  of  Equation  (64)  it  was  seen  that  aj  =  ( 1 / A ) f (aa)ln(B/A)  where 
B/A  is  defined  in  Equation  (63).  This  means  that  Eoi (l/A)ln(B/A)  =  1  or 
bj  =  f ( ct3 ) / A  becomes 


bi 


f(aa) 


EqiI  n 


EEqi  -  Eq2\ 

Eq2  ‘  Eqi  J 


(87) 
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In  Equations  (86)  ana  (87)  me  values  cf  Eq2  depenas  on  the  magnitude  of 
me  applied  load  as  it  must  ae  evaluated  after  the  removal  of  the  first 
I oaa  appl i cat i on  . 

Altnougn  this  section  defines  the  dependence  of  the  parameters  of 
Euqations  1 6  7 )  and  (54)  on  the  load,  it  may  be  more  expedient  to 
determine  6t  ana  m  for  Ecuation  (67)  or  31  ana  pj  for  Equation  (64)  from 

experimental  ooservations  on  repeated  loadings.  in  this  case,  the 
experiment  must  be  made  using  me  vertical  and  confining  stress  expected 
in  the  field.  The  vertical  stress  is  predictable  from  stress  transfer 
theory  as  shown  in  Equation  (76).  and  the  confining  pressure  is  the 
product  of  the  coefficient  of  earth  pressure  at  rest,  the  effective  unit 
weight,  ana  me  depth. 
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SECTION  VI 

RUTTING  MODEL  FOR  MULTILAYERED  PAVEMENTS  WITH  PARTICULATE  MATERIAL 

A.  INTRODUCTION 

The  designs  against  rut: .  r,g  are  currently  either  empirical  or 
quasi -el asti c.  The  empirical  method  selects  pavement  thicknesses  based 
cn  correlations  of  excessive  deformations  to  subgraae  strength  or  to 
satisfy  an  allow sole  value  of  the  vertical  subgrade  strain.  The 
subgrade  strength  is  usually  taken  as  its  California  Bearing  Ratio  (CBR) 
value  or  more  recently  recommended  is  its  resilient  modulus  (Thompson, 
1984).  These  methods  cannot  be  used  to  predict  the  amount  of 
deformation  present  after  load  applications.  The  quasi  -  el asti c  method 
is  more  direct  and  has  tne  capability  of  oDtaining  cumulative 
deformations.  This  approach  uses  elastic  theory,  either  linear  or 
nonlinear,  to  predict  the  expected  stress  state  within  the  pavement. 
The  approach  is  termed  "quasi -el asti c”  since  it  uses  this  predicted 
elastic  stress  state  to  design  laboratory  tests  to  measure  the 
relationship  of  permanent  strains  to  number  of  repeated  load 
appl i cations . 

In  linear  elastic  analysis  the  pavement  is  assumed  to  be  composed 
of  homogeneous,  isotropic  linear  elastic  layers  that  are  infinite  in  the 
horizontal  extent.  The  loading  is  vertical  and  uniformly  distributed  on 
a  circular  base  with  no  surface  shear.  It  is  also  assumed  that  there  is 
full  continuity  between  the  layers  and  within  the  layers.  The  main 
limitation  of  this  approach  is  that  the  aggregate  base  of  the  pavement 
is  not  linear  elastic.  In  nonlinear  elastic  analysis,  the  same 
assumptions  are  made  except  that  the  elastic  modulus  of  the  granular 
layer  is  described  as  a  function  of  the  stress  state  based  on  an 
empirical  curve  that  fits  experimental  results.  This  approach  requires 
complex  numerical  modeling,  and,  through  its  requirement  for  full 
continuity  between  the  layers  and  among  the  aggregate  and  voids, 
violates  the  no  tensile  stress  requirement  of  the  aggregate. 
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Designs  based  primarily  on  the  stress  transfer  of  layered  linear 
or  nonlinear  elastic  theory  fail  to  evaluate  the  benefits  of  properly 
graded  and  compacted  granular  material.  In  general,  the  method  of 
compaction  and  gradation  of  the  granular  material  affects  the  nature  of 
tne  stress  distribution  in  the  material  and  if  properly  controlled  can 
be  used  to  increase  the  carrying  capacity  of  the  pavement  by  spreading 
the  load,  thereby  reducing  stresses  in  the  lower  layers  and  decreasing 
the  stress  concentration  that  causes  rutting.  It  was  already  shown  in 
Section  II  that  the  particulate  theory  of  stress  transfer  includes  a 
diffusion  coefficient  that  is  a  function  of  the  method  of  compaction  and 
gradation  of  the  granular  material.  It  is  shown  in  this  section  how 
this  method  of  stress  transfer  can  oe  modified  to  predict  stresses  in 
layered  pavements  with  particulate  and  elastic  layers.  The  strain 
accumulation  theory  of  Section  III  is  also  extended  to  multiple  layers 
in  order  to  predict  rutting  in  the  pavement. 

B.  STRESS  TRANSFER  IN  MULTILAYERED  MEDIA 


Flexible  pavements  generally  consist  of  a  2-  to  8-inch  asphalt 
surface  over  a  12-  to  36-inch  granular  base  or  subgrade.  Therefore,  it 
is  a  three-layer  system,  with  each  layer  transmitting  stress  in  a 
different  manner  due  to  the  difference  in  the  type  of  material.  In  this 
section,  the  particulate  stress  transfer  theory  of  Section  II  will  be 
modified  to  represent  the  condition  of  stress  transfer  through  pavement 
layers  of  varying  properties. 

The  general  expression  for  the  vertical  stress  induced  by  an 
inclined  point  load  of  vertical  component  0  and  horizontal  component  qiQ 
in  the  X  direction  is  given  by  Equation  (31)  and  can  be  rewritten  as 


q  f-  p2  -2qtx  +  qt2  _  i 

°z  4nW ( z  )  exp  1  4W( z )  J 


(88) 


where  p2  =  x2  +  y2. 
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To  determine  the  vertical  stress  under  an  inclined  uniform  load  q 
on  a  circular  area  of  radius  "a",  Equation  (88)  can  be  modified  by 
replacing  0  Dy  qrdraO,  p  d>  r2  +  p2  -  2rpcos0  ana  x  by  x  -  r  cos  (£+0) 

where  £,  =  tan_I(y/x).  The  final  equation  is 


az  = 


Q 


4rcw  ( 2 ) 


-2n  /’a  |r2+p2-2rp  cos0-2a[X+2qjrcos(0+^)+q[222 

exp! - r4wT7j - ; - 


rdrde  (89) 


For  the  special  case  of  a  vertical  load  (q^O),  this  equation  can 
be  integrated  over  r  to  give 
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COS20/W( Z  ) 


For  known  values  of  "a"  and  W(z).  this  latter  equation  can  be 
easily  integrated  numerically  at  any  value  of  p.  Directly  under  the 
center  of  the  vertically  loaded  area  (p=0),  Equation  (90)  gives 


az  =  q 


exp 


-  a2 
_4W(z) 


(91) 


Equations  (89),  (90),  and  (91)  present  progressively  simpler 
cases.  However,  in  each  of  these  equations,  all  of  the  geometric  and 
material  properties  are  contained  in  the  function  W(z).  Therefore,  the 
determination  of  the  stress  distribution  in  layered  media  requires  only 
a  modification  of  W(z)  in  that  layer.  Directly  below  the  center  of  a 
vertically-loaded  circular  area  of  radius  "a"  on  granular  material  it 
was  shown  in  Section  IV  that  W(z)  =  Wg(z),  where  Wg(z)  refers  to 

granular  layer  and  is  given  by 

Wg(z)  =  di  z2  +  d2az  (92) 
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-sere  aj  anu  g->  are  materi al - sceci f i c  constants  reflecting  the  particle 
sizes,  oacKir'.c.  ana  previous  '.jading.  An  equivalent  function,  W(z),  can 
also  oe  found  for  material s  that  transmit  stress  close  to  that  predicted 
Dy  elastic  theory.  In  this  case,  setting  Equation  (77)  equal  to 
Equation  (75)  gives  W(z)  =  We(z),  wrtere 


We(z) 


(93) 


In  order  to  determine  tne  modification  of  W(z)  in  a  multilayered 
pavement,  tne  method  presented  by  Golden  (1984)  is  used.  In  this 
method,  W-j(z)  represents  the  function  W(z)  for  stress  distribution  in  a 

media  composed  only  of  the  material  in  the  ith  layer,  and  the  function 
W ( z )  in  Equation  (89),  (90).  or  (91)  is  taken  as 

W i ( z )  for  z  1  hi 

W(z)  =  -  W2(z)  -  W2(hi)  +  W i ( h i )  for  hi  <  z  1  hi  +  h2  (94) 

W3 ( z )  •  W3 ( h i-*-h2 )  +  W 2 ( h  1  h 2 )  ■  W2 ( h  1 )  +  W 1  ( h  1 )  for  z  2.  hj+h2 

to 

where  h i i s  the  thickness  of  the  ith  layer. 

The  stress  distribution  in  the  asphalt  layer  is  elastic  in  nature. 
Therefore,  the  value  of  W 1 ( z )  for  flexible  pavements  is  given  by 

Equation  (93).  The  second  layer  of  the  pavement  is  the  granular  layer, 
and  the  value  of  W2 ( z )  is  given  by  Equation  (90).  The  third  layer  is 

the  subgrade.  For  cohesive  subgrade  it  was  observed  that  due  to  its 
cohesive  nature  the  elastic  prediction  of  stresses  is  acceptable  (Morgan 
and  Gerrard,  1981).  In  this  case.  Equation  (93)  gives  W3 ( z ) .  For  sand 

subgrade,  the  elastic  solution  is  not  acceptable  (Morgan  and  Gerrard, 
1981),  and  W3 ( z )  is  of  the  form  of  Equation  (92).  The  vertical  stress 
on  the  subgrade  (z  =  hi  +  (12)  directly  below  the  center  of  the  loaded 
area  is  found  from  Equation  (91),  where  from  Equation  (94) 
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W(z)  =  di  (  t\%  +  2 h | h2 )  +  d?  ah2  +  - 


6  In 


a 

[TU  1 


+  1 


(95) 


The  role  playea  oy  me  layer  thicknesses  and  material  constants  on 
distributing  the  stress  to  the  weaker  subgrade  layer  is  apparent  from 
tnis  equation.  From  Equation  (89)  a  larger  W(z)  translates  into  smaller 
vertical  stresses.  Therefore,  increases  in  all  quantities  h j .  h2.  di . 
and  d2  decrease  the  stress  transferred  to  the  weaker  subgrade.  However, 
since  parameters  di  and  d2  are  functions  of  the  gradation  and  compaction 
of  the  granular  layer.  Equation  (95)  can  be  used  to  evaluate  the 
benefits  of  increased  layer  thicknesses  versus  increased  compaction  and 
gradation . 

C.  RUTTING  PREDICTION  MODEL 

The  majority  of  research  on  subgrade  and  granular  materials  has 
shown  a  log-log  relationship  of  the  permanent  strain  to  the  number  of 
load  applications  (Yoder  and  Witczak,  1974).  This  log-log  relationship 
was  also  derived  theoretically  for  asphalt  materials  (Khedr,  1988)  and 
observed  experimental ly  (Majidzadeh.  Khedr,  and  El  -Moj arrish.  1979). 
This  relationship  is  identical  to  Equation  (67)  derived  for  granular 
media  with  m  not  equal  to  one. 

For  granular  material,  it  was  also  observed  that  for  values  of  m 
approaching  one.  Equation  (64)  becomes  the  more  appropriate  model.  The 
parameter,  Bi  in  Equation  (67)  (and  ai  in  Equation  (64))  is  the 
permanent  deformation  after  the  first  cycle  of  loading  and  for  granular 
material  is  given  by  Equation  (54).  In  this  sense,  Biis  an  explicit 
function  of  the  vertical  stress  at  the  point  in  question  in  the  media. 
The  parameter  m  is  also  a  function  of  the  stress  level,  as  shown  in 
Section  V.  Therefore,  it  is  imperative  that  in  the  evaluation  of  Bi  and 

m  that  the  correct  stress  be  used. 
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The  tneoretical  result  of  Section  III  ana  those  of  Khedr  (1986) 
indicate  t.nat  at  any  depth  z  under  the  loaded  area  the  accumulated 


permanent  strain  after  N  cycles  is  given  py 


P 

£(  z ) 
aN 


.■ilc?) 

B(<*z)N 


(96) 


The  parameters  B(az)  ana  m(az )  are  functions  of  the  vertical 
stress  ana  material  characteristics.  If  subscripts  a,  g,  and  s  refer  to 
the  asphalt  layer,  granular  layer,  ana  subgrade,  then  the  permanent 
depression  in  the  flexible  pavement  after  h  cycles  is 


5n 


= 

(’hi  ma(oz) 

3a (az )  N  dz+ 

fhl+ho  ma(cz) 

3g(Oz)  N  y  dz+ 

J 

0 

hi 

fTlc  ( O7  ) 

Bs(<*z)  N  dz  (97) 


hi+h2 


Equation  (97)  is  the  general  prediction  model.  Although  it  is  possible 
to  express  the  B(az)  parameters  as  continuous  known  functions  of  z.  this 
is  not  possible  for  the  m(az)  parameters.  This  is  illustrated  in 
Section  V  for  the  granular  layer  and  in  Khedr  (1986)  for  the  asphalt 
layer.  As  a  result,  the  integrals  of  Equation  (97)  must  be  replaced  by 
summations,  and  the  pavement  divided  into  discrete  sections  of  magnitude 
Az .  The  values  of  B(ctz)  and  m(cz)  are  then  evaluated  using  the  stress 

expected  in  the  center  of  each  section.  Further,  although  it  is  shown 
in  Section  V  that  ’t  i:  poss'ole  to  express  B(az)  and  m(crz)  in  terms  of 

other  material  parameters,  it  is  more  efficient  to  determine  these 
parameters  directly  from  repeated  load  tests  using  linear  regression  on 

p 

the  log-log  transformation  of  the  measured  values  of  ean  and  N.  This 
method  concentrates  on  the  evaluation  of  only  two  parameters  B  and  m  per 
laboratory  test.  However,  it  is  imperative  that  the  stress  state 
expected  in  the  field  be  used  to  perform  the  laboratory  test.  As  an 
improvement  over  the  use  of  stresses  determined  from  elastic  theory,  the 
values  of  B(crz)  and  m(az)  should  be  evaluated  using  the  stress  predicted 

by  the  multilayered  theory  described  above.  A  simple  closed-form 
description  of  the  vertical  stress  to  be  applied  to  the  laboratory 


samel e  is  provided  by  the  ccmDinaiion  of  Equations  (91)  and  (94)  where 
i-i  nz)  and  W3  ( z )  are  aescrioea  tv  Equation  (91)  and  W2  ( z )  by 

Equation  (92). 
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SECTION  VII 

CONCLUSIONS  AND  RECOMMENDATIONS 


CONCLUSIONS 


The  report's  objective  is  to  present  the  results  of  research 
investigating  rutting  in  flexible  pavements  from  a  perspective  whereby 
the  granular  layer  is  treated  as  particulate  rather  than  a  continuum. 
It  snows  that  tne  stress  transfer  in  the  granular  layer  can  be  derived 
without  tne  assumption  of  the  existence  of  the  second  derivation  of 
strains.  This  assumption  is  generally  used  for  both  linear  and 
nonlinear  elastic  stress  analysis.  However,  in  the  granular  layer,  the 
particles  are  discrete  and  deformation  is  tne  result  of  discrete 
particle  motion  rather  than  continuous  particle  compression. 

From  the  results  of  the  report,  the  following  conclusions  can  be 

drawn : 

•  A  particulate  definition  of  stress  and  strain  in  granular 

media  is  necessary  in  order  to  adequately  model  the  behavior  of  the 
granular  layer  in  flexible  pavements.  This  definition  is  necessary 
since  stresses  are  transferred  only  at  particle  contacts  and  not  through 
the  voids  in  the  material,  and  strains  are  the  result  of  particle 

movements  rather  than  particle  compression. 

•  The  particulate  definition  of  stress  and  strain  allows  the 

development  of  a  stress  continuity  relationship  which  when  combined  with 
the  equilibrium  equation  gives  the  stresses  in  the  material.  This 
development  has  three  major  advantages  over  conventional  linear  and 
nonlinear  continuum  analyses.  First,  the  assumption  that  the  second 
derivatives  of  strains  exist  at  all  points  is  not  unrealistically 
imposed  on  the  discrete  medium;  secondly,  the  observed  stress-strain 
behavior  of  the  material  is  built  into  the  stress  continuity 

relationship;  and  thirdly,  the  stresses  are  expressed  in  terms  of  a 

material  term  W(z)  that  is  a  function  of  the  gradation  and  method  of 

compaction  of  the  medium. 
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•  r-.e  : :  r:  ’  cu  ;  a  te  theory  prediction  of  stress  and  strain 

conforms  v  exceri mental  ooservations  and  ail  material  constants  can 
ze  evaluated  i  tr.  tne  use  of  conventional  experimental  methods. 

•  'r:e  particulate  theory  of  stress  distribution  reduces  to  the 

elastic  case  ..’tr  the  appropriate  choice  of  W(z),  and  the  evaluation  of 
stress  transfer  tr.rougn  multilayered  media  is  easily  determined  with  a 
systematic  representati on  of  W(z)  in  terms  of  the  thicknesses  and 
materials  composing  the  layers. 

•  The  particulate  approach  to  modeling  granular  behavior 

confirms  tne  log-log  relationship  between  accumulated  permanent  strain 
and  number  of  cycles  of  loading.  However,  it  magnifies  the  fact  that 
parameters  ir,  this  relationship  are  highly  dependent  on  the  stress  state 
at  tne  point  wnere  the  strain  is  evaluated. 

•  The  particulate  theory  of  stress  distribution  through 

multilayered  media  can  be  used  to  predict  the  stress  state  to  be  used 
for  evaluation  of  the  material  constants  in  the  permanent  strain 

accumulation  versus  number  of  cycles  of  loading.  This  is  an  improvement 
over  conventional  linear  and  nonlinear  elastic  analysis  due  to  its 
simplicity  and  the  rationality  of  the  particulate  approach  to  stress 
,  distribution  in  granular  media. 

I 


B.  RECOMMENDATIONS 

The  evaluation  of  rutting  in  flexible  pavements  presented  in  this 
report  is  based  on  a  new,  more  rational  definition  of  stress  and  strain 
in  particulate  media.  This  definition  assumes  that  deformation  is  the 
result  of  particle  movement  rather  than  particle  compression  and  that 
stress  is  transferred  only  at  particle  contacts.  It  is  gratifying  that 
this  definition  very  simply  leads  to  models  that  adequately  represent 
observed  s tres s - s t ra i n  behavior  for  granular  soils  under  one 
dimensional,  isotropic,  triaxial.  and  shear  loading  conditions  (see 
Appendix  A).  It  also  easily  illustrates  the  diffusive  nature  of  stress 
transfer  in  granular  material  under  inclined  loading.  Therefore,  it  is 


57 


■  e commenced  zr.ii  c  trier  •  col ems  concerr~j  „un  the  transmission  of 
:  tress  ::!g  strain  r.  granule  v  m  e  d  i  a  o-„-  aoproached  -with  this  new 
:-ir ;  •  i :  1  on  .  Included  among  tr.ese  are  tr.e  fc'i  ’owing: 

•  Modeling  tne  cerformance  of  sana  grids  under  loads.  Sand 

a d  s  are  ..sea  increasingly  for  expeditionary  airfields  and  their  design 
.,  i '  1  ce  greatly  improved  if  the  mecnani  sm  of  tneir  behavior  is 
•'at'onally  modeled.  This  approach  is  ideal  for  such  investigation. 

•  Developing  improved  methods  of  interpreting  the  information 

obtained  from  nondestructive  testing  of  r'lexiole  pavement  using  devices 
such  as  tne  falling  weignt  cef 1 ectometer .  This  new  approacn  will  allow 
,ror  tne  evaluation  of  more  characteristic  pavement  properties  rather 
tnan  the  resilient  modulus  currently  measured. 

•  investigating  the  mechanics  of  the  nonlinear  nature  of 

stress  wave  transmission  in  granular  soils,  thereby  improving  the 
understanding  of  Dlast  ana  earthquake  effects  on  structures  and 
potentially  liquefiable  soil. 
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A  GENERAL  STRESS-STRAIN  MODEL  FOR  GRANULAR  SOILS 


R i n g s i e v  Harrop- Williams 


The  SDM  Corp.,  McLean,  VA  22102 


\3STR ACT 

'>  general  stress-strain  model  is 
ienvea  for  granular  soils  based  on  the 
premise  that  the  deformation  in  these 
soils  is  not  caused  so  much  by  the 
ie formation  of  individual  particles  as  by 
the  relative  movements  of  the  more  mobile 
particles.  The  model  is  general  and 
specific  relationsnips  are  obtained  for 
different  loading  conditions.  For  one 
dimensional  and  isotropic  loading 
conditions  a  stress-strain  curve  concave 
towards  the  stress  axis  is  derived,  and 
for  triaxial  and  shear  loading  conditions 
the  derived  s t r e ss- s t r a i n  curve  is  concave 
towards  the  strain  axis.  In  particular 
soils  that  show  a  distinct  peak  stress  are 
well  modeled,  and  the  model  reduces  to  the 
hyperbolic  model  for  soils  exhibiting  no 
distinct  peak  stress.  All  cases  show 
excellent  fit  to  experimental  data. 


INTRODUCTION 

The  stress-strain  behavior  of  granular 
soils  is  created  primarily  by  individual 
particle  movements  to  form  a  denser  matrix 
rather  than  elastic  compression  of  the 
particles.  As  a  result  the  experimental 
curve  is  always  nonlinear. In  general  it  is 
concave  towards  the  strain  axis  under 
triaxial  and  shear  loading,  and  concave 
towards  the  stress  axis  in  one  dimensional 
and  isotropic  loading  conditions.  To 
adequately  predict  the  settlement  under  a 
load  one  needs  to  be  able  to  model  this 
nonlinear  behavior. 

Of  primary  importance  in  settlement 
calculations  is  the  need  to  model  the 
triaxial  compression  condition.  One  simple 
empirical  method  proposed  by  Konder  and 
Zelasko  (4)  is  the  widely  used  hyperbolic 
model.  This  model  appears  to  be  a  natural 
fit  for  soils  that  do  not  exhibit  a 
distinct  peak  stress.  Another  empirical 
model  proposed  is  a  modification  of  the 
Ramberg-Osgood  model  used  for  dynamic 
loading  for  static  conditions  (2).  This  is 
based  on  fitting  a  curve  to  the  initial 
tangent  mjuulus,  the  modulus  of  the 
plastic  zone,  the  yield  stress  and  a 


parameter  defining  the  order  of  the  curve. 
For  order  one  this  curve  becomes  a 
hyperbola.  One  other  approacn  suggested  is 
to  treat  the  tangent  shear  modulus  and 
tangent  bulk  modulus  as  variables  (5). 
These  are  usually  taken  as  linear 
functions  of  the  octahedral  normal  and 
shear  stresses. 

The  difficulty  in  modelling  t.ne 
nonlinear  stress-strain  behavior  of  sand 
under  load  has  led  many  investigators  to 
propose  numerical  curve  fitting 
techniques.  One  of  the  most  popular  of 
these  is  the  piecewise  linear  method.  Here 
the  nonlinear  experimental  curve  is 
divided  into  pieces  of  linear  elastic 
sections  for  numerical  analysis.  Very 
often  these  are  the  incremental  Hooke's 
law  or  the  hypoelastic  law  (2).  Another 
numerical  method  is  the  use  of  spline 
functions  to  fit  experimentally  observed 
curves.  These  are  functions  that  uses  the 
data  to  provide  an  analytic  curve  similar 
to  the  graphical  process  of  using  a  French 
curve  (1).  These  also  require  intense 
numerical  procedures,  and  the  data  must  be 
presented  in  a  smooth  form  and  not 
scattered  as  observed  experimentally. 

In  this  paper  a  general  stress-strain 
model  is  derived  from  basic  principles. 
This  model  is  shown  to  apply  to  one 
dimensional,  isotropic,  triaxial,  and 
shear  loading  conditions  as  special  cases 
based  on  the  boundary  conditions  they 
impose. 

STRESS  AND  STRAIN  IN  GRANULAR  SOILS 

In  the  determine  cion  of  stresses  and 
strains  in  granular  soils  any 
representative  element  of  the  soil  with 
volume  dxdvdz  must  be  composed  of 
particles  and  voids.  The  element  must  also 
consist  of  enough  particles  such  that  the 
relative  movements  of  the  particles  as  a 
result  of  forces  on  these  particles 
produce  strains.  This  means  that  in  the 
limit  the  volume  dxdvdz  cannot  be  made  to 
approach  zero  but  must  instead  approach 
some  minimal  finite  volume  ijh.  This 
minimal  volume  of  fixed  dimensions  i,  j 
and  h  in  the  x,  y  and  z  directions, 
respectively,  is  the  smallest  volume  of 
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granular  soil  that  can  as  used  to  define 
-cress  ana  strain.  In  effect,  cnis  volume 
.3  an aiojous  to  a  point  in  continuous 
•catena!  .  The  latter  r. avina  a  rixea  volume 
:  f  cero  magnitude  on  a  macroscopic  scale 
cut  it  i  microscopic  level  is  composea  of 
:  iscrete  atoms  that  enter  ana  leave  cr.e 
point. 

Since  particles  are  allowea  to  enter 
and  leave  the  minimal  element,  if  w  is  the 
average  relative  displacement  of  a 
particle  in  the  element  in  the  z 
iirection,  then  the  normal  strain  in  the 
element  i.n  the  z  direction  is 

£  .  =  w  /  h  1 ) 


direction.  Also  let  k  r  represent  the 
i  v  e  r  a  a  e  resistance  of  a  particle  to 
movement  in  tie  x  direction.  Therefore, 
the  force  in  the  x  direction  in  the  lower 
half  of  element  1  is  F,  i  =  k,  u  ,  and  the 
force  in  the  x  direction  in  the  upper  half 
of  element  2  is  rj2  =  N';  kz  [  u  +  (  3u/ ?z  )  dz  ]  . 
Here  S;  is  the  numoer  of  particles  in  the 
lower  half  of  element  1  and  S';  is  the 
number  of  particles  in  the  upper  half  of 
element  2.  We  assume  that  S  does  not 
change  much  with  z  as  compared  with 
changes  in  u,  therefore,  N i  =  N;  =  N/2. 

The  snear  stress  in  the  x  direction  at  the 
interface  between  the  two  elements  is 
(Ftl-rx2)/(dxdv),  or  since  dxdvdz  approaches 
ijh  in  the  limit  the  snear  stress  is 


For  simoiicitly  we  introduce  a  microscopic 
stiffness  coefficient  k ,  to  represent  the 
average  resistance  of  a  particle  to 
movement  in  the  z  direction  sucn  that  the 
force  on  a  particle  in  the  minimal  element 
: n  the  z  direction  is  ktv.  The  magnitude 
of  k j  depends  on  the  rougnness  of  the 
particles  and  the  confining  pressure.  The 
force  in  the  element  in  the  z  direction  is 
Fz  =  Nkz w ,  where  N  is  the  number  of 
particles  in  the  element.  The  normal 
stress  in  the  element  in  the  z  direction 
is  then  Fz / (  i  j) 

Oj.  =  N'kz  w/(  i  j)  (2) 

In  this  equation  N'  =  ijh/[Vp(i  +  e)],  where 
Vo  is  the  average  volume  of  a  particle  and 
e  is  the  void  ratio.  By  definition  the 
volumetric  strain  is  related  to  the  void 
ratio  as  e. =  (  e, -e ) / ( 1+e,  )  ,  where  e0  is 
the  initial  void  ratio.  Therefore, 

V  =  ijh/[Vp ( 1+e, )(!-£„)]  (3) 

The  substitution  of  equations  fl)  and  (3) 
into  equation  (2)  and  taking  the  z 
direction  as  the  axial  direction  give  the 
general  relationship  of  axial  stress  to 
axial  strain  as 

3i-  E0 e  a  /  ( ( 1  -  S £»)  (4) 

where 

to  =  h*  kz  /  [  Vp  (  1+e,  )  ]  (  5  ) 

and  the  parameter  3  was  added  since  the 
condition  »  1  is  unattainable  under 
conventional  loads.  From  equation  14),  the 
derivative  of  Ca  with  respect  to  £«  at 
ev  =  ea=  0  is  E0.  Therefore,  E,  is  the 
initial  tangent  modulus  of  the  soil. 

To  look  at  shear  stresses  in  granular 
soils  consider  two  minimal  elements  tf 
soil  adjacent  to  each  other  but  separated 
by  a  surface  dxdy  in  the  xv  plane.  Let 
element  1  be  at  location  z  with  average 
particle  displacement  u  in  the  x 
direction,  and  element  2  be  at  location 
z+dz  with  average  relative  particle 
displacement  u+13u/3z)dz  in  the  x 


r  =  [ Nk* h / ( 2 i j ) ] (3u/3z  )  (6) 

ONE  DIMENSIONAL  AND  ISOTROPIC  COMPRESSION 


In  one  dimensional  compression  of  soils 
e»  =  £s  ,  and  in  isotropic  compression  we 
have  £v =  3£a  .  From  equation  (4)  these  two 
conditions  can  be  represented  by  the 
single  equation 

oa=  E,£a  / ( 1-  nea)  (7) 

where  0  is  the  reciprocal  of  the 
assvmptotic  strain  observed  at  large 
stress.  Figure  (1)  shows  the  excellent  fit 
of  equation  (7)  to  data  for  isotropic 
loading  conditions  on  two  samples  of 
McCormic  Ranch  sand  obtained  from 
reference  (2,  pg.  193). 
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Fig.  1.  Stress-Strain  in  Isotropic  Loading 
(1  psi  =  6.895  KPa ) 

One  should  also  note  that  equation  (7) 
can  be  written  as  aa  »  E0  (  Ga  +  HEJ+  .  .  ) 

and  by  series  expansion  we  have  that 
fexpi2n£a)-i]/(2h)  =  £a  +  nej+  .  .  . 
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There  tore,  letting  a  =  2  n  a  n  a  A  =  E  0 / a , 

: .ves  the  alternate  expression 

a„  =  Aexp(a  ea )-A  { S) 

Equation  (  8  )  is  the  same  as  that  derived 
by  the  nvpoeiastic  analysis  tor  these 
loading  conditions  (2,  pg  139). 

TRIAXIAL  COMPRESSION 

In  triaxial  compression  of  elastic 
material  the  radial  strain, e.,  is 
proportional  to  the  axial  strain,  and  the 
proportional  constant  is  the  empirical 
Poisson's  ratio.  In  effect,  if  v.  =  -  d £r ,/d£. 
and  vs =  -  £r/£a  ,  then  in  elastic  material 
v. •  Vs  .  However,  in  triaxial  compression 
of  sands  the  relationship  of  radial  strain 
to  axial  strain  is  nonlinear  (3).  To 
represent  this  nonlinear  condition  one  can 
let  v, -  v.  =  D,  where  D  is  an  empirical 
constant  representing  the  average 
difference  between  J,  and  v,.  If  D  =  0  the 
relationsnip  is  linear,  and  the 
nonlinearity  increases  with  the  magnitude 
of  D.  The  definition  v,  =  _  £r /  £a  gives 
d  v, /d£,  =  (tr  +  £av,)/tj  =  —  D/£» ,  also 
since  the  volumetric  strain  £«  =  £,  +  2£r 
we  have  vs  =  (  £a-£v)/(2£.  i,  which  shows 
that  dv,/d£,  =  -0.5d(s. /  £a )/de, .  Equating 
the  expressions  for  dv,/d£a  gives  the 
expression  d(G*/£«)/dG.  =  2D/G»  ,  which  has 
solution 

£»  =  BE.  +  2D  e.lne.  (9) 

where  B  is  an  integration  constant.  It 
should  be  noted  that  xlnx  approaches  zero 
as  x  approaches  zero. 

The  relationship  of  volumetric  strain 
to  axial  strain  under  triaxial  compression 
for  a  medium  dense  sand  at  two  different 
confining  pressures  as  obtained  from 
reference  (2,  pg.  177)  is  shown  in  Figure 
(2a).  The  regression  of  equation  (9)  on 
the  data  in  thi  '  figure  is  shown  as  the 
dashed  lines.  An  excellent  fit  is  observed 
for  each  case. 

The  substitution  of  equation  (9)  into 
equation  (  4  )  yields 

Oa  =  Eo  £*/fl-3BC*  -23D£aln£.)  (10) 

Letting  aa  ”  and  £a  -  £p  at  maximum  stress, 
we  get  6B  «  1 /e? -  Eo/Op  -  28DlnEp.  Also, 
setting  dO»/d£,  =  0  at  maximum  stress 

gives  2BD  =  -1/Cp  .  The  substitution  of 
these  into  equation  (10)  gives  the 
relationship 

0,  =  £,  /  [  a  f  (  £.  )  +  be.  ]  (11) 

where  a  =  1/E„,  b  =  l/Op,  and 

fee.)  =  i  4-  ( ca  /  gp)  [  ln(  £a  /£p  )  -  1]  (12) 
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pressure  £p  approaches  infinity.  In  this 
case  equation  (12)  gives  f(e»)‘.  1,  and 
equation  (11)  reduces  to  the  hyperbolic 
model 

Oa=£a/(a+bea)  (13) 

The  excellent  fit  of  equation  (11)  to  the 
s t r e s s - s t r a i n  responses  of  a  medium  dense 
sand  at  two  different  confining  pressures 
is  shown  in  Figure  (2b). 

SHEAR  LOADING 

The  shear  strain  in  the  element  of 
soil  during  shear  loading  is  Y=  3u/9z. 
Therefore,  the  substitution  of  equation 
(3)  into  equation  (6)  gives 


In  soils  with  no  distinct  peak  stress  such 
as  loose  sands  and  sands  under  high 


t  -  Go  Y/(  1-  6ev) 


(14) 


where 

Go  =  h*  kx  /  [  2V?  (  I  +  e0  J  1  (15) 

ina  the  parameter  3  was  added  since  the 
condition  £v  =  l  is  unattainable  under 
conventional  loads.  Evaluating  dT/dY  at 
v  =  e,  =  o  shows  that  G0  is  the  initial 
tangent  shear  modulus. 

In  elastic  material  under  shear  loading 
the  volumetric  strain  is  proportional  to 
the  shear  strain.  This  means  that  in 
•lastic  material  d£v/dY  =  £„/Y .  In 
granular  soils,  however,  this  linear 
relation  does  not  hold  (6).  To  represent 
the  nonlinear  relationship  between  ev  and Y 
we  let  i£„,'dY  -  £v/Y  =  C,  where  C  is  an 
empirical  constant  representing  the 
average  difference  between  the  values 
:£v/JY  and  £v/Y  .  This  expression  says  that 
£„  =  v [ ( d£v /  dY  ) -C  ]  .  Differentiating  both 
sides  of  this  with  respect  to  Y  gives 
i-£„,dY“  =  C/Y,  wnich  has  solution 

£v  =  B,Y  +  cYlnY  (16) 

vnere  Bn  is  an  integration  constant. 

Substituting  equation  (16)  into 
equation  ( 14)  and  letting  T  =  Tp  and  Y  =  Yp 
at  peak  stress  gives  the  first  constant  as 
B Bn  —  1,'Yp  -  Go/fp  -  SClnYp.  Further, 
setting  df/dY  =  0  at  peak  stress  gives  the 
constant  BC  =  -  1/Yp  .  The  substitution  of 
these  into  equation  (14)  gives  the 
expression 

t  =  Y / [ af (Y  )  +  bY]  (17) 

where  a  =  1/Go,  b  =  1/Tp,  and 

f(Y)  =  1  a-  (Y/Yp  )[ln(Y/YP)  -  1]  (18) 

“or  soils  with  no  distinct  peak  stress  we 
have  the  condition  Yp  =  «•  and  equation  (18) 
gives  f(Y)  =  1.  In  this  case  equation  (17) 
reduces  to  the  hyperbolic  model. 

The  direct  shear  experimental  data  for 
the  change  in  sample  thickness,  t,  with 
respect  to  shearing  displacement,  s,  for  a 
sand  in  a  loose  and  a  dense  state  as 
presented  by  Sutton  (6,  pg.  99),  are  shown 
in  Figure  (3a).  The  corresponding  shear 
stress  versus  displacement  data  for  this 
sand  is  shown  in  Figure  (3b).  The  sand  was 
tested  under  drained  conditions  and  a 
constant  normal  stress  of  150  kN/m2  was 
applied  by  the  piston  in  both  tests.  In 
the  direct  shear  test  the  change  in  sample 
thickness  t  is  proportional  to  the 
volumetric  strain  £, and  the  shearing 
displacement  s  is  proportional  to  the 
shear  strain  Y  .  Therefore,  the 
relationship  of  t  to  s  is  of  the  same  form 
as  equation  (16).  That  is 

t  =  Bj  s  +  sins  (19) 

where  3;  and  Cj  are  constants.  Also,  the 
relationsnio  of  T  to  s  is  of  the  same  form 
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(b)  Shear  Stress  Versus  Displacement 
Fig.  3.  Direct  Shear  loading 
as  equations  (17)  and  (18),  or 

T=s/[a1f(s)+b^s]  (20) 

and 

f(s)  =  1  +  ( s/sp  )  [ ln( s/sp )  -  1]  (21) 

where  sp  is  the  displacement  at  peak 
stress,  and  and  b[  are  constants. 

The  regression  of  equation  (19)  on  the 
data  of  Figure  (3a)  are  shown  as  the 
dashed  linos  of  the  figure.  Also  the 
regression  of  equation  (20)  on  the  data  of 
Figure  (3b)  are  shown  as  the  dashed  lines 
on  that  figure.  In  each  case  the  fit  is 
excellent . 


SUMMARY 


Based  on  the  assumption  that  the 
stress-strain  behavior  of  granular  media 
is  controlled  by  the  displacements  of 
individual  particles  rather  than  particle 
compression  a  general  stress-strain  model 
for  granular  material  is  derived.  This 
model  says  that  for  both  axial  and  shear 
loading  the  stress  is  proportional  to 
r.  s  t  r  a  in  /  (  l-6ev)  ]  ,  where  the  proportional 
constant  is  the  initial  modulus  of  the 
soil  and  3  is  a  constant.  For  one 
dimensional  and  isotropic  loading 
conditions  the  volumetric  strain  £v  is 
proportional  to  the  axial  strain.  For 
triaxial  and  shear  loading  conditions 
nonlinear  relationships  of  e« to  the  axial 
strain  and  £*  to  the  shear  strain  are 
developed.  It  is  shown  that  the  model  can 
be  used  for  soils  that  show  distinct  peak 
stresses  when  under  triaxial  and  shear 
loading,  and  that  it  reduces  to  the 
hyperbolic  model  for  soils  that  shows  no 
distinct  peak  stress. 
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